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An  experimental  investigation  was  conducted  into  the 

effects  of  surface  roughness  on  the  losses  in  the  flow 

j  ■  ■  ■  ■  •  d  •’  .  . 

through  a  two-dimensional  compressor  cascade^*  The  cascade 
consisted  of/  seven  NACA  65-A506  airfoils  with  two  inch  , 

chords  and  aspect  ratios  of,  one ^  The  blades  were  mounted 

with  a  stagger  angle  of  16  degrees,  an  angle  of  attack  of  15 

/ 

degrees  and  a  solidity  of  1.5.  The  Reynolds  number  pel:  foot 
was  in  excess  of ,»fcwo  million.  Three  blade  configurations 
were  evaluated:  one  smooth  and  two  with  different  degrees 

of  roughness  applied  to  the  first  quarter  chord.  Hot 
wire/film  anemometers  were  used  -to*  measure*/ velocities  and 
turbulence  intensities  both  in  the  wake  and  along  the  suc¬ 
tion  surface  of  the  center  blade.  Pressure  distributions 
were  obtained  from  taps  drilled  along  the  suction  surface. 

Surface  roughness  is  shown  to  Increase  both  boundary 
layer  edge  velocities  and  boundary  layer  thickness  towards 
the  trailing  edge  while  advancing  the  transition  point 
towards  the  leading  edge.  It  is  also  shown  to  Increase  wake 
full  thickness,  momentum  thickness,  form  factor,  and  total 
pressure  loss  coefficient  while  decreasing  the  velocity  re¬ 
covery  immediately  downstream  of  the  trailing  edge.  The 
wake  characteristics  are  most  sensitive  to  an  initial  in¬ 
crease  in  roughness,  but> regardless  of  the  roughness  level, 
the  majority  of  the  changes  occur  within  the  first  half  *a 
chord  length  distance  downstream. 
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SURFACE  ROUGHNESS:  ITS  EFFECTS  ON  THE 


PERFORMANCE  OF  A  TWO-DIMENSIONAL 
COMPRESSOR  CASCADE 

I .  Introduction 

The  nature  o£  the  £low  through  a  rotating  compressor 
blade  row  Is  a  complicated  one.  Factors  Including  viscos¬ 
ity,  unsteadiness,  and  three-dimensionality  make  it  dif¬ 
ficult  to  obtain  purely  theoretical  solutions.  Therefore,  in 
order  to  obtain  the  necessary  blade  design  data,  the 
"...prevailing  approach  has  been  to  treat  the  flow  across 
individual  compressor  blade  sections  as  a  two-dimensional 
flow"  (13:183).  This  lends  itself  to  experimental  investi¬ 
gations  using  two-dimensional  cascades.  Though  some  recent 
studies  have  questioned  the  reliability  of  using  cascade 
data  in  predicting  turbomachine  blade  performance  (8:292), 
it  is  "...nevertheless  felt  that  the  calculation  and  under¬ 
standing  of  the  flow  through  two-dimensional  cascades  will 
remain  the  corner-stone  of  turbomachlnery-aerodynamics  for  a 
long  time  to  come"  (24:496). 

The  Air  Force  Institute  of  Technology  (AFIT),  in  con¬ 
junction  with  the  Air  Force  Wright  Aeronautical  Labora¬ 
tories,  has,  over  the  past  few  years,  made  use  of  a  cascade 


to  study  the  effects  of  surface  roughness  on  the  performance 
of  a  compressor.  Interest  In  this  area  is  based  on  the 
knowledge  that  compressor  efficiency  is  influenced  a  great 
deal  by  the  surface  quality  of  the  blading  (3:283).  A  num¬ 
ber  of  examples  exist  in  the  literature  (Refs  3  and  22) 
where  full  scale  testing  of  compressors  were  used  to  collect 
data.  Acknowledging  that  this  is  the  preferred  method  of 
analysing  component  performance,  it  is  still  felt  that  exa¬ 
mination  of  the  flow  through  a  two-dimensional  cascade  can 
contribute  to  this  study. 

Research  began  in  1981  with  the  design  and  construction 
of  the  AFIT  cascade  facility.  Since  then,  successive  inves¬ 
tigations  (Refs  17,25,30)  have  improved  the  facility,  exa¬ 
mined  the  variation  of  total  pressure  loss  with  roughness, 
and  initiated  surveys  of  the  cascade  blade  pressure  profile 
and  boundary  layer. 

Objectives  and  Scope 

The  objectives  of  this  Investigation  were:  (1)  to  con¬ 
centrate  bn  the  suction  surface  of  the  center  blade  in  the 
cascade  and  obtain  more  detailed  information  on  the  effects 
of  surface  roughness  on  its  pressure  profile  and  boundary 
layer  growth;  (2)  examine  the  flow  in  the  center  blade's 
wake  to  determine  trends  in  the  variation  of  certain  wake 
parameters  with  roughness;  and  (3)  quantify  the  losses  due 
to  roughness  experienced  by  the  cascade  in  terms  of  a  total 
pressure  loss  coefficient. 
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Hot  wire/£ilm  anemometers  were  used  to  measure  velo¬ 
cities  and  turbulence  intensities,  and  the  center  blade 
suction  surface  was  instrumented  with  pressure  taps  to  ob¬ 
tain  the  pressure  profile.  A  secondary  goal  of  this  study 
was  to  fully  automate  the  control  and  movement  of  the  probe 
around  the  suction  surface  contour  and  integrate  that  fea¬ 
ture  into  the  data  acquisition  system.  Surface  roughness 
was  measured  by  a  profilometer  and  characterized  in  terms  of 
an  arithmetic  average  roughness,  Ra  ,  an  equivalent  sand 
roughness,  ks  ,  and  a  so-called  technical  roughness  k  .  The 
latter  was  formulated  by  Schaffler  (22:10)  in  an  attempt  to 
better  quantify  the  roughness  produced  by  different  kinds  of 
manufacturing  processes.  It  is  presented  as  a  comparison 
and  not  used  in  describing  the  results  of  this  study. 
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Evaluating  the  losses  o£  a  flow  through  a  cascade  re¬ 


quires  a  knowledge  of  the  boundary  layer  around  the  blades, 
the  pressure  distribution  about  the  contour  of  the  blades, 
and  the  effects  of  mixing  in  the  wake  of  the  cascade 
(23:771) . 

Boundary  Laver 

"Although ...  the  blade  surface  boundary  layer  constitutes 
a  small  portion  of  the  flow  field,  it  plays  a  decisive  part 
in  the  determination  of  the  actual  flow  characteristics  of 
the  cascade"  (21:153).  This  is  particularly  true  in  the 
case  of  the  suction  surface  of  a  blade  where  the  presence  of 
an  adverse  pressure  gradient  is  conducive  to  flow  separa¬ 
tion. 

when  analyzing  boundary  layer  velocity  profiles,  one 
must  take  into  account  the  effects  of  the  normal  pressure 
gradient  that  exists  between  the  suction  surface  of  one 
blade  and  the  pressure  surface  of  the  next.  With  the  static 
pressure  remaining  nearly  constant  across  a  boundary  layer 
(21:152),  Bernoulli's  equation  yields  a  velocity  at  the  edge 
of  the  pressure  surface  boundary  layer  that  is  less  than 
that  at  the  corresponding  suction  surface  location.  Mea¬ 
suring  the  velocities  between  the  two  blades  would  therefore 
yield  a  profile  similar  to  that  depicted  in  Fig.  1(a).  As  a 
result  of  this  effect,  the  boundary  layer  edge  velocity  and 


thickness  are  difficult  to  determine.  One  way  to  approach 
the  problem  is  to  assume  that  the  measured  profile  is  actu¬ 
ally  a  composite  profile  consisting  of  three  regions.  These 
are:  (1)  an  inner  region  close  to  the  surface  where  vis¬ 
cosity  acts;  (2)  an  outer  invlscid  region  in  the  channel 
where  the  normal  pressure  gradient  acts;  and  (3)  a  region 
between  the  two  where  the  conditions  match  (7:8).  This  can 
be  represented  mathematically  through  the  use  of  "The  Method 
of  Matched  Asymptotic  Expansions"  detailed  by  Van  Dyke 
(28:77-97).  Essentially,  the  matching  principle  used  is 
that  the  inner  limit  of  the  outer  region  is  equal  to  the 
outer  limit  of  the  inner  region.  The  composite  profile  can 
then  be  obtained  by  summing  the  inner  and  outer  expansions 
minus  the  part  they  have  in  common  so  as  not  to  count  it 
twice.  Mellor  and  Wood  (16:307),  in  their  study  of  axial 
compressor  end-wall  boundary  layer  theory,  expressed  this  as 
follows : 

Une*a  =  Uw.  +  Uihv  ~  Ut  (  1  ) 

where  Umeas  *3  the  measured  velocity  of  the  composite  pro¬ 
file,  the  boundary  layer  velocity,  the  velocity 
in  the  invlscid  region,  and  Ug  the  boundary  layer  edge 
velocity.  The  three  regions  of  the  composite  profile  are 
shown  in  Fig.  1(b).  As  a  result  of  this  procedure,  the  edge 
velocity  can  be  determined  as  well  as  the  boundary  layer 


thickness . 


Pressure  Distribution 

The  variation  In  local  static  pressure  around  the  con¬ 
tour  of  a  blade  is  usually  presented  as  a  dimensionless 
pressure  coefficient  Cp  .  It  is  defined  as 

C,  =  (PL  -  Pi )  /  ( 1/2  PV,1)  (2) 

where  PL  is  the  local  static  pressure,  and  P1  ,  P  ,  am 
Vj  the  static  pressure,  density,  and  velocity,  respec¬ 
tively,  at  the  inlet  to  the  cascade. 

Information  obtained  from  the  pressure  distribution 
include  the  detection  of  flow  separation,  observed  as  a 
flattening  of  the  pressure  profile.  In  certain  cases,  a 
separated  laminar  boundary  layer  may  become  turbulent  and 
reattach  itself  to  the  surface.  This  forms  a  region  of 
recirculating  flow  known  as  a  laminar  separation  bubble 
(23:220).  Roberts,  as  reported  by  Scholz  (24:475-476), 
described  the  formation  of  these  bubbles  by  first  distin¬ 
guishing  between  two  regions  in  the  flow.  The  first  lies 
between  the  separation  streamline  and  the  surface  and 
actually  makes  up  the  bubble.  The  other  is  the  free  shear 
layer  itself  extending  from  the  same  streamline  out  to  the 
inviscid  flow.  While  this  shear  layer  remains  laminar,  the 
flow  is  separated  and  represented  by  a  constant  pressure 
plateau  in  the  pressure  profile.  If  the  shear  layer  becomes 
turbulent,  the  entrainment  of  high  energy  free-stream  fluid 
allows  the  boundary  layer  to  reattach  itself  to  the  surface. 
Mueller  (18:3),  in  his  investigation  of  separated  flow  near 


the  leading  edge  of  airfoils.  Interpreted  sharp  discontin¬ 
uities  In  the  pressure  profile  following  the  plateau  as  the 
location  of  the  onset  of  transition  and  the  point  of  reat¬ 
tachment  . 

Flow  In  the  Wake 

The  boundary  layers  that  form  on  the  suction  and  pres¬ 
sure  surfaces  of  a  blade  in  a  cascade  meet  at  the  trailing 
edge  and  form  a  wake.  As  this  wake  moves  downstream,  the 
retarded  flow  is  mixed  with  the  free-stream  until,  at  some 
distance  behind  the  blade  trailing  edge,  the  flow  Is  once 
again  uniform  (21:153).  This  Is  pictured  In  Fig.  2. 

The  velocity  defect  In  the  wake  also  manifests  itself  as 
a  local  defect  in  total  pressure.  The  mixing  that  occurs 
downstream  represents  an  additional  loss  of  energy  so  that 
the  total  pressure  In  the  plane  of  uniform  flow  Is  less  than 
that  at  the  trailing  edge  (15:4). 

Total  Pressure  Loss  Coefficient .  Though  the  effects  of 
the  cascade  on  the  flow  can  be  expressed  in  many  ways,  in¬ 
cluding  the  various  force  coefficients  such  as  drag,  it  is 
best  described  in  terms  of  the  loss  in  total  pressure.  In 
this  form,  it  is  "...especially  adapted  to  the  practical 
demands  of  turbomachine  design"  (24:34).  This  loss  para¬ 
meter,  known  as  the  total  pressure  loss  coefficient  3  ,  is 
defined  to  be  the  ratio  of  the  mass-averaged  loss  in  total 
pressure  across  the  blade  row  to  a  reference  free-stream 
dynamic  pressure  (13:201).  In  this  investigation,  it  was 


obtained  from  the  following: 
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where  P,p  is  the  total  pressure  upstream  of  the  cascade,  P 

the  mass-averaged  total  pressure  across  a  blade  spacing 

2 

downstream,  and  1/2  P  the  dynamic  head  at  the  inlet. 

Unlike  some  parameters  such  as  the  drag  coefficient,  13 


will  vary  with  distance  downstream  of  the  trailing  edge  due 
to  the  effects  of  viscosity  and  turbulent  mixing.  Conse¬ 
quently,  the  complete  loss  in  total  pressure  through  the 
cascade  is  one  based  on  the  properties  of  the  flow  at  a  sta¬ 
tion  far  enough  downstream  so  that  the  flow  is  once  again 
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uniform  across  the  blade  spacing  (21:156). 

Wake  Thicknesses  and  Parameters .  As  can  be  seen  from 
Fig.  2,  the  flow  field  for  some  distance  downstream  of  the 
trailing  edge  of  a  blade  can  be  separated  into  a  wake  region 
and  a  free-stream  region.  It  would  therefore  be  of  interest 
to  express  the  total  pressure  los3  in  terms  of  the  wake 
characteristics  (15:5). 

Like  a  boundary  layer  over  a  surface,  various  wake 
thicknesses  can  be  calculated  from  the  velocity  distribution 
normal  to  the  outlet  flow.  The  actual  wake  thickness  <5  is 
arbitrarily  defined  to  be  the  width  of  the  wake  as  esta¬ 
blished  by  the  points  where  V/Vo  are  equal  to  0.99 
(14:15).  Making  reference  to  Fig.  3,  the  displacement 
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thickness  <3  ^  ,  momentum  thickness  62  •  and  pseudoenergy 
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The  velocity  profile  In  the  wake  can  also  be  character¬ 
ized  by  so-called  shape  factors.  These  are  dimensionless 
quantities  formed  by  taking  the  ratio  of  the  thicknesses 
discussed  above  (23:673).  Two  of  particular  use  in  analyz¬ 
ing  the  wake  are  the  form  factor  H  and  the  pseudoenergy 
factor  K  defined  as  follows  (15:10): 


<3i/<5. 

(7) 

<5*/ <3, 

(8) 

Wake  Minimum  Velocity.  The  downstream  variation  of  the 
wake  minimum  velocity  is  one  way  of  measuring  the  intensity 
of  the  velocity  defect  in  the  wake.  Lieblein  and  Roudebush 
(14:6-7)  postulate  that  this  variation  is  in  the  form 

V„„,/V0  -  1  -  a(x/c  +  b)',/*  (9) 

where  vmin  1s  the  minlmum  velocity  in  the  wake,  V0  the 
free-stream  velocity,  x  the  distance  downstream  from  the 
trailing  edge,  and  c  the  chordlength  of  the  cascade  air¬ 
foil.  They  found  the  value  of  the  constants  a  and  b  to 
be  0.13  and  0.025,  respectively.  They  arrived  at  this  con¬ 
clusion  by  comparing  actual  turbulent-wake  cascade  data  to 
theoretical  studies  done  on  laminar-wake  profiles  behind 
solid  bodies. 

Wake  Form  Factor .  Lieblein  and  Roudebush  (14:8)  also 
determined  that  an  empirical  variation  of  form  factor  with 
downstream  distance,  obtained  previously  for  isolated  air¬ 
foils,  is  also  applicable  for  the  turbulent  wakes  of  cascade 


airfoils.  The  relationship  Is 
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where  Htb  is  the  form  factor  at  the  trailing  edge  of  the 

blade  and  b  the  same  empirical  constant  found  in  Eq  (9) 

Wake  Full  Thickness .  The  downstream  variation  of  the 
wake  full  thickness  ratio,  <5/c  ,  for  isolated  airfoils,  as 
derived  from  empirical  data,  is  given  by 


(3/0,  =  (<J/c)r«  +  0 . 052  ( x/c  )e-n 


(11) 


where  (3  /c)<j*gls  the  full  thickness  ratio  at  the  blade  trail¬ 
ing  edge.  Lieblein  and  Roudebush  (14:11-12)  believe  that 
this  relationship  al3o  holds  for  cascades  due  to  the  simi¬ 
larities  observed  between  the  mixing  processes  in  the  wakes 
of  cascade  and  isolated  airfoils. 

Total  Pressure  Loss .  A  theoretical  analysis  of  total 
pressure  loss  in  low-speed  two-dimensional  cascade  flow  was 
carried  out  by  Lieblein  and  Roudebush  (15).  In  it,  they 
formulate  the  following  equation  for  the  total  pressure  loss 
coefficient  at  an  arbitrary  distance  x  downstream  of  the 
trailing  edge: 


■"  11  *  9-H> 


(12) 


where  and  Oy^  are  the  angles  between  the  flow  direc¬ 

tion  and  the  cascade  axial  direction  at  the  inlet  and  par¬ 
ticular  outlet  plane,  respectively.  The  wake  momentum 
thickness  parameter  8^  included  in  the  equation  is  given 
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(3»/c)i  cr/cos  13% 


<13) 


where  or  is  the  solidity  of  the  cascade. 

In  order  to  obtain  the  theoretical  loss  equation  in  its 
presented  form,  an  approximate  relationship  between  the 
pseudoenergy  factor  K  and  the  form  factor  H  was  used. 
Values  for  K  were  obtained  for  different  analytical  repre¬ 
sentations  of  wake  velocity,  and  it  was  determined  that  the 
K-H  relationship  based  on  a  power  velocity  profile  was  both 
simple  and  accurate  (15:11-12).  It  was 
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K  =*  (H  +  1 ) /  ( 3H  -  1) 


(14) 


As  discussed  previously,  the  complete  loss  in  total 
pressure  is  measured  at  a  station  far  enough  downstream  so 
that  the  velocity  Is  uniform  across  the  blade  spacing.  The 
total  pressure  loss  coefficient  at  this  station,  Om  ,  can 
be  expressed  in  terms  of  the  wake  characteristics  at  an  ar¬ 
bitrary  measuring  plane  downstream  of  the  trailing  edge  by 


cos_01l*  1  +  1/2  6,H,» 
COS  /3i  J  "  (1  -  e,H,)* 


(15) 


This  equation,  developed  by  Lieblein  and  Roudebush  (15:22), 
Incorporates  the  analytical  observation  that  for  values  of 
less  than  1.2,  the  air  outlet  angle  has  very  little 

effect  on  the  loss  for  complete  mixing. 
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Cascade  Test  Facility 

The  experimental  work  for  this  thesis  was  conducted  on 
the  Cascade  Test  Facility  (CTF)  located  in  the  Air  Force 
Institute  of  Technology  School  of  Engineering.  The  CTF  con¬ 
sisted  of  a  40  horsepower  centrifugal  blower  with  a  dis¬ 
charge  rating  of  3000  ft/'min  at  26  oz  of  head.  Outside  air 
was  dtawn  into  a  sealed  containment  room  where  the  option 
existed  to  mix  it  with  recirculated  warmer  room  air  in  order 
to  stabilize  the  discharge  temperature.  Upstream  of  the 
blower  itself  were  screen  wire  filters  and  an  electrostatic 
air  cleaner  assembly  for  trapping  particles.  The  outlet  flow 
from  the  blower  was  directed  through  a  nine  foot  long  dif¬ 
fuser  section  into  a  stilling  tank  where  it  was  radially 
diffused  by  a  center  body  plug.  Before  exiting  the  tank,  the 
flow  was  filtered  once  more  by  a  40  mesh  wire  screen  and 
then  straightened  by  a  four  inch  thick  honeycomb  grid.  It 
then  entered  a  two  by  eight  inch  test  section  with  a  Rey¬ 
nolds  Number  per  foot  in  excess  of  two  million  and  inlet 
turbulence  intensities  generally  less  than  two  percent  as 
measured  by  a  hot  film  sensor. 

A  more  detailed  description  of  the  CTF  is  given  by 
Allison  (1).  A  general  schematic  of  the  hardware  is  shown 


INLET  DUCT  TEST  SECTION 


Test  Section 

The  CTF  test  section,  shown  in  Fig.  5,  consisted  of 
seven  NACA-A506  airfoils.  According  to  Moe  (17:5),  the  pro¬ 
file  of  these  blades  was  similar  to  that  found  in  the  latter 
stages  of  a  high  pressure  compressor.  The  blades  were  set 
at  a  row  angle  of  31  degrees,  a  stagger  angle  of  16  degrees, 
and  an  angle  of  attack  of  15  degrees.  This  resulted  in  a 
flow  turning  angle  of  19  degrees.  The  individual  blades  had 
a  chord  and  span  of  two  inches  for  an  aspect  ratio  of  one. 
The  solidity  was  1.5  based  on  a  blade  spacing  of  1.333  in. 
In  addition,  the  two  outside  blades  were  half  imbedded  into 
the  end  walls  in  order  for  the  test  section  to  simulate  an 
infinite  cascade  (30:8). 

Seven  static  pressure  taps  were  located  on  one  side  of 
the  inlet  duct  to  the  test  section  and  were  used  to  verify 
the  uniformity  of  the  upstream  flow  and  determine  its  velo¬ 
city.  The  exit  channel  was  made  up  of  two  13  in  long  plexi¬ 
glass  sides  with  adjustable  endwalls.  Four  static  pressure 
taps  were  distributed  along  the  length  of  each  endwall  to 
aid  in  setting  the  channel  to  the  desired  configuration. 
Additional  pressure  taps  were  placed  along  one  sidewall  to 
coincide  with  the  four  exit  plane  traverse  locations. 

S ldewa 1 1  Suction  System 

As  described  by  Moe  (17:8),  a  boundary  layer  control 
mechanism  was  incorporated  into  the  CTF  test  section.  The 
sidewalls  immediately  around  the  blades  were  made  of  a 


porous  stainless  steel  through  which  the  boundary  layer  can 
be  continuously  drawn  off.  Suction  was  provided  by  an  In¬ 
dustrial  vacuum  cleaner  connected  to  the  test  section  by  two 
aluminium  manifolds.  Moe  (17:47)  determined  that  with  the 
suction  applied,  two  dimensional  flow  was  established  over 
the  center  span  of  the  airfoils  in  the  cascade  extending  for 
approximately  two  thirds  the  width  of  the  blade. 

Instrumentation 

Monitoring  of  the  flow  conditions  in  the  CTF  required 
the  use  of  a  variety  of  pressure,  temperature,  and  velocity 
measuring  devices.  A  complete  listing  of  the  components 
that  make  up  these  systems  is  attached  as  Appendix  A.  Mea¬ 
surements  made  include  the  stilling  tank  total  pressure  and 
temperature,  test  section  inlet  and  exit  static  pressures, 
and  ambient  temperature  and  pressure.  Eighteen  30-ln  U-tube 
water  manometers  were  also  integrated  into  the  system  to 
monitor  the  test  section  inlet,  exit,  and  endwall  pressure 
taps  as  well  as  the  tank  total  pressure  and  sidewall  suc¬ 
tion.  In  addition,  a  48  channel  Scanivalve  pressure  measur¬ 
ing  system  was  used  to  obtain  the  pressure  distribution 
about  the  suction  surface  of  the  center  blade  of  the 
cascade  . 

The  pressure  transducers  were  calibrated  over  their 
operating  ranges  with  a  linear  curve  fit  used  to  convert  the 
output  voltages  to  pressures.  The  two  "T-type"  copper/con- 
stantan  thermocouple  voltage  outputs  were  converted  into 


temperatures  through  the  use  of  software  provided  with  the 
data  acquisition  system  (10:8.78,8.85). 

Velocity  and  turbulence  data  in  the  blade  boundary  layer 
and  wake  regions  were  obtained  using  three  Thermal  Systems 
International  (TSI)  Model  1050  Constant  Temperature  Anemo¬ 
meters.  Two  were  used  with  a  TSI  model  1241-10  "X"  config¬ 
uration  hot  film  probe  to  survey  the  center  blade's  wake  re¬ 
gion.  The  remaining  one  was  connected  to  a  TSI  model 
1218-T1.5  hot  wire  boundary  layer  probe.  Since  the  temper¬ 
ature  of  the  air  going  through  the  CTF  is  increased  by  30  to 
40  o  F  above  ambient,  the  hot  wire/film  probes  were  calibra¬ 
ted  using  a  scheme  that  takes  into  account  the  effects  of 
elevated  temperatures  (29).  The  complete  procedure  is  de¬ 
tailed  in  Appendix  B.  The  final  result  is  a  single  calibra¬ 
tion  curve  for  each  sensor  that  is  applicable  over  the  range 
of  expected  test  conditions.  Figures  6  and  7  are  examples 
of  such  curves  for  the  wake  and  boundary  layer  probes,  re¬ 
spectively. 

The  calibration  procedure  required  the  use  of  an  air 
supply  whose  velocity  and  temperature  could  be  measured  and 
varied  as  needed.  This  was  accomplished  using  a  modified 
TSI  Model  1125  Calibrator  system  instrumented  to  measure  the 
total  temperatures  at  the  interior  chamber  inlet  and  exit  as 
well  as  the  total  pressure.  A  portion  of  the  source  air, 
provided  by  the  shop  compressor,  was  diverted  through  a 
quartz  heater  assembly  powered  by  an  adjustable  20  amp 
transformer.  This  "hot"  air  was  then  mixed  with  "cool"  air 
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directly  from  the  compressor  before  being  directed  into  the 

calibrator.  By  varying  the  amount  of  "hot"  and  "cool"  air 
mixed  upstream,  the  velocity  and  temperature  of  the  air 
exiting  the  calibrator  could  be  adjusted.  A  specified  tem¬ 
perature  was  reached  when  the  two  thermocouples  agreed  to 
within  a  tenth  of  a  degree  Fahrenheit. 


Traversing  System 

The  hot  wire/film  sensors  used  to  probe  the  flow  about 
the  cascade  were  connected  to  an  18  in  probe  support,  the 
opposite  end  of  which  was  fitted  in  an  airfoil  shaped  mount. 
The  complete  sensor/support  assembly  was  attached  to  a  tra¬ 
versing  mechanism  capable  of  motor  driven  movements  in  both 
the  X  and  Y  directions.  Spanwise  movement  was  possible  by 
manually  adjusting  a  thumbscrew.  The  traversing  system  as  a 
whole  could  be  rotated  about  the  Z  axis  in  order  to  align 
the  probe  to  the  desired  angle  in  relation  to  the  cascade. 

Though  the  sensor  could  be  positioned  manually,  the  tra¬ 
versing  mechanism  was  integrated  into  the  data  acquisition 
system  to  allow  automatic  control  of  probe  movement.  In  this 
mode,  the  accuracy  of  sensor  positioning  was  estimated  at 
+/-  0.002  in  in  both  the  X  and  Y  directions.  The  actual  re¬ 
solution  of  the  stepper  motors  in  the  X  and  Y  directions 
were  0.002  and  0.001  in,  respectively. 


Blade  Roughness  Configurations 

The  NACA  65-A506  airfoils  used  in  this  investigation 
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were  cast  ln  sets  of  six  blades  each  using  Fiber-Re3in  Cor¬ 
poration  FR44  resin  and  5595  hardener.  The  mixture  was 
poured  into  an  aluminium  casting  mould  and  cured  at  room 
temperature  for  ten  hours.  To  further  increase  the  stiff¬ 
ness  of  the  blades,  the  castings  were  then  heat  treated  in 
two  hour  increments  at  150,  200,  and  250  °F. 

It  can  be  deduced  that  roughness  would  have  the  greatest 
effect  on  the  performance  of  a  cascade  if  It  were  applied  on 
the  suction  surface,  near  the  leading  edge  of  the  blade.  In 
this  region,  the  boundary  layer  is  subjected  to  an  adverse 
pressure  gradient  and  the  relative  magnitude  of  the  rough¬ 
ness  elements  to  the  boundary  layer  thickness  is  largest. 
Tanis  (25:50)  demonstrated  this  experimentally  by  varying 
the  location  of  surface  roughness  and  measuring  the  effect 
on  the  pressure  loss  coefficient. 

In  previous  investigations  using  the  AFIT  CTF  (17;25;30) 
the  primary  method  of  adding  roughness  to  a  blade  was  by 
applying  different  grades  of  carborundum  grit  to  the  surface 
and  sealing  them  ln  place  with  an  acrylic  adhesive.  This 
procedure  had  a  number  of  drawbacks,  the  most  serious  of 
which  was  the  possibility  of  altering  the  blade  profile, 
particularly  when  applying  the  larger-sized  grit.  To  allev¬ 
iate  this  problem,  the  application  of  surface  roughness  in 
this  study  was  incorporated  into  the  actual  casting  process. 
This  was  done  by  tacking  a  strip  of  emery  paper  to  the  suc¬ 
tion  side  of  the  casting  mould,  beginning  a  sixteenth  of  an 
inch  back  from  the  leading  edge  and  extending  to  the  25  per- 


cent  chord  point.  Two  grades  of  emery  paper  were  used:  280 
and  180.  In  addition,  one  set  of  blades  were  cast  without 
added  roughness  to  serve  as  a  baseline  smooth  configuration 
for  comparison. 

The  surface  roughness  of  each  blade  was  measured  with  a 
Rank  Taylor  Hobson  Surtronic  3  prof i lometer .  A  cut-off  dis¬ 
tance  of  0.8  mm  was  used  with  the  standard  pickup  112/1503. 

A  total  of  36  readings  were  taken  about  the  center  blade's 
midspan  and  averaged.  The  output  of  the  profilometer  was 
the  arithmetic  average  roughness,  Ra  ,  in  micrometers.  This 
roughness  parameter  Is  described  in  Appendix  C  along  with 
two  empirical  correlations  to  convert  Ra  to  the  equivalent 
sand  roughness,  ka  ,  and  Schaffler's  technical  roughness, 
k  .  The  results  are  tabulated  below. 


TABLE  I 

Blade  Roughness  Data 


Con£.  # 

Ra,  4m 

ka, 

k,  4* 

1 

0.45 

2.79 

4.01 

2 

12.1 

75.0 

107.7 

3 

18.3 

113.5 

162.9 

For  comparison,  Koch  and  Smith  (12:423-424)  report  an  Ra 
value  of  18  micrometers  for  180  grade  sandpaper.  For  the 
remainder  of  this  investigation,  roughness  will  be  charac¬ 
terized  in  terms  of  Ra  or  k«j  . 
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Suction  Sur face  Pressure  Distribution 

Information  on  the  suction  surface  pressure  profile  was 
obtained  by  instrumenting  the  center  blade  of  each  roughness 
configuration  with  38  static  pressure  taps.  Due  to  the 
small  size  of  the  blade,  the  taps  were  offset  a  sixteenth  of 
an  inch  on  either  side  of  the  midspan  and  alternated  left 
and  right.  The  use  of  a  staggered  layout  for  the  taps  was 
possible  due  to  the  two-dimensionality  of  the  flow.  In  ad¬ 
dition,  Williams  (30:24)  demonstrated  that  a  profile  ob¬ 
tained  with  a  staggered  arrangement  did  not  differ  signifi¬ 
cantly  from  one  obtained  with  all  the  taps  aligned  along  the 
blade's  centerline. 

Holes  were  drilled  from  either  end  of  the  blade  to  join 
with  the  taps.  Small  (0.022  in  OD )  metal  tubes  were  inser¬ 
ted  into  the  holes  and  connected  to  the  Scanivalve  system 
with  pressure  lines.  The  actual  pressure  tap  arrangement  is 
depicted  in  Fig.  8.  A  complete  profile  was  not  possible  due 
to  the  thin  trailing  edge,  and  the  need  to  accommodate  holes 
for  the  blade's  four  mounting  pins. 

Data  Acquisition  and  Anal vs i3  System 

The  CTF  was  monitored  and  controlled  from  an  air  condi¬ 
tioned  room  located  alongside  which  housed  a  Hewlett-Packard 
(HP)  9845B  computer  integrated  with  an  HP  3052A  Automatic 
Data  Acquisition  System.  Again  the  various  components  of 


the  system  are  listed  in  Appendix  A. 

Software  controlled  the  movement  of  the  traversing  mech- 


anism  and  read  the  voltage  outputs  from  the  various  measur¬ 
ing  devices.  Though,  in  certain  cases,  data  were  reduced 
on-line,  the  standard  procedure  used  was  to  store  the  vol¬ 
tage  readings  on  magnetic  disks.  These  were  later  converted 
and  stored  as  engineering  units  making  them  available  to  ad¬ 
ditional  programs  used  to  calculate  the  performance  parame¬ 


ters  of  interest. 


IV.  Experimental  Procedure  and  Data  Reduction 

Testing  Procedure 

The  testing  procedure  for  this  study  was  divided  Into 
three  stages.  The  first  stage  consisted  of  a  survey  of  the 
center  blade's  wake  with  a  1241-10  hot  film  probe.  This  was 
followed  by  the  measurement  of  the  pressure  distribution 
about  the  suction  surface  of  the  center  blade.  The  third 
and  final  stage  made  use  of  a  1218-T1.5  boundary  layer  probe 
to  survey  the  flow  field  around  the  suction  surface  of  that 
same  blade.  This  three  part  procedure  was  repeated  for  each 
configuration  of  roughness. 

Prior  to  taking  any  measurements,  the  CTF  was  allowed  to 
warm  up  and  the  flow  to  reach  a  stabilized  operating  temper¬ 
ature.  The  exit  endwalls  and  the  boundary  layer  suction 
were  then  adjusted  to  balance  the  test  section.  The  ba¬ 
lanced  condition  was  reached  when  the  static  pressure  taps 
along  both  endwalls  read  approximately  ambient  pressure,  and 
the  inlet  pressure  to  the  test  section  was  uniform.  Moe 
(17)  demonstrated  that  balancing  the  test  section  in  this 
manner  established  two  dimensional  flow  through  the  cascade. 

Wake  Survey 


With  the  test  section  balanced,  the  traversing  mechanism 
was  rotated  about  the  Z  axis  so  as  to  align  the  probe  sup¬ 
port  with  the  exit  flow  direction.  This  direction  was  de- 


termined  from  the  position  of  the  endwalls.  wake  surveys 
were  then  conducted  at  1.25,  2.25,  3.25,  and  4.25  in  down¬ 
stream  of  the  center  blade's  trailing  edge.  Each  traverse 
consisted  of  133  data  points  located  at  0.01  in  intervals  in 
the  Y  direction  as  shown  in  Fig.  9. 

Velocity  Correction .  The  voltage  outputs  of  the  two 
sensors  were  converted  into  X  and  Y  components  of  velocity 
through  the  use  of  the  probe's  calibration  curve.  As  was 
found  in  previous  work  with  this  cascade  facility 
( 25  : 16; 17 : 12; 30  :  23  ) ,  the  measured  velocity  was  five  to  ten 
percent  higher  than  theoretically  possible.  This  was  ap¬ 
parent  when  the  exit  total  pressure  calculated  from  the 
measured  velocity  was  higher  than  the  pressure  measured  in 
the  stilling  tank.  This  discrepancy  was  thought  to  be  re¬ 
lated  to  the  variation  in  heat  transfer  rate  with  varying 
humidity  as  well  as  the  variation  of  the  probe  support  re¬ 
sistance  with  changing  flow  temperature  (25:14-15).  Note 
that  during  the  calibration  procedure,  the  probe  support  is 
not  exposed  to  the  different  temperatures  of  the  flow  to  the 
extent  that  it  is  during  an  actual  run.  The  problem  was 
corrected  using  a  procedure  outlined  by  Tanis  (25:15-16)  and 
Moe  (17:12)  in  which  continuity  between  two  centerline 
planes  located  upstream  and  downstream  of  the  cascade  is 
used . 

Local  Turbulence  Intensity.  The  ratio  of  the  root  mean 
square  velocity  to  the  streamwise  velocity  at  a  particular 
data  point  was  used  as  the  definition  of  the  local  turbu- 
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TU  =  Vmm/V,  (16) 

where  V  is  the  local  time  varying  velocity,  measured  as 

KfIS 

the  RMS  voltage  output  of  the  anemometer,  and  is  the 

local  axial  velocity,  measured  as  the  DC  output. 

Due  to  the  configuration  of  the  two  sensors  on  the 
probe,  Tu  in  the  wake  could  be  resolved  into  X  and  Y  compo¬ 
nents,  thereby  giving  an  orientation  to  the  fluctuations. 

Total  Pressure  Loss .  The  total  amount  of  losses 
through  a  cascade  was  quantified  by  the  nondimensional  total 
pressure  loss  coefficient,  defined  earlier  in  Eq  (3).  The 
total  pressure  upstream  of  the  cascade  was  measured  by  the 
pressure  transducer  situated  in  the  stilling  tank.  The  to¬ 
tal  pressure  at  a  particular  traverse  location  downstream 
was  calculated  by  first  determining  the  total  pressure  P  at 
each  data  point  as  follows: 


(V,»  )* 

*2c7fT  J 


(17) 


where 


P,j  *  exit  static  pressure 
Vt'  “  corrected  exit  velocity 
Tt  *  tank  total  temperature 
y  *  ratio  of  specific  heats 
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The  exit  total  pressure  could  then  be  determined  by  mass 
averaging  these  individually  calculated  total  pressures: 
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where  PT^  is  the  mass-averaged  exit  total  pressure,  p  the 
local  density  evaluated  using  the  perfect  gas  law,  and  dA 
the  incremental  area. 

As  pointed  out  by  Moe  (17:16),  the  area  integrals  above 
can  be  reduced  to  single  Integrals  with  respect  to  Y  since 
the  flow  is  two-dimensional  and  spanwise  uniform  at  the 
blade  centerline.  The  integrals  were  solved  numerically 
using  the  trapezoidal  rule. 


Suction  Surface  Pressure  Distribution 

The  Scanivalve  pressure  measurement  system  wa3  used  to 
measure  the  static  pressure  at  each  of  the  38  tap  locations 
on  the  suction  surface  of  the  center  blade.  This  was  re¬ 
peated  up  to  five  times  for  each  run,  and  the  results  stored 
on  magnetic  disks.  The  pressure  coefficients  were  calcula¬ 
ted  using  Eq  ( 2 ) . 

This  measurement  was  carried  out  Immediately  prior  to 
any  boundary  layer  run.  This  allowed  for  the  calculation  of 
a  local  velocity,  at  each  tap  location  based  on  the 
local  static  pressure.  Applying  Bernoulli's  equation  along 
a  streamline  from  the  inlet  of  the  test  section  to  a  partic¬ 
ular  tap  location  yields 


This  local  velocity  was  used  as  a  check  against  the 
results  obtained  with  the  boundary  layer  probe  during  the 
next  phase  of  the  investigation.  In  addition,  it  was  used 
to  calculate  the  mean  reference  temperature  at  which  the 
fluid  properties  would  be  evaluated  (see  Appendix  B). 


Suction  Surface  Boundary  Layer  Survey 

The  boundary  layer  along  the  midspan  of  the  center 
blade's  suction  surface  was  surveyed  at  thirteen  chord  lo¬ 
cations.  These  corresponded  to  pressure  taps  and  were  lo¬ 
cated  at  the  4.69,  9.38,  25,  29.69,  34.38,  40.63,  45.31,  50, 
65.63,  70.31,  75,  79.69,  and  84.38  percent  chord  positions. 

In  the  previous  investigation  (30:17),  limitations  in 
the  axial  movement  of  the  traversing  mechanism  prohibited 
measurements  forward  of  the  42.75  percent  chord  location. 
This  limitation  was  resolved  in  this  study  by  modifying  the 
probe  support  mount  fairing  to  allow  an  additional  2.25  in 
of  movement  in  the  axial  direction.  This  modification  also 
alleviated  the  blockage  problem  experienced  by  Williams 
(30:17),  negating  the  requirement  to  rebalance  the  test  sec¬ 
tion  for  each  traverse. 

Probe  Positioning  and  Control .  One  of  the  objectives  of 
this  thesis  was  to  fully  automate  the  acquisition  of  boun¬ 
dary  layer  data.  This  required  that  the  movement  of  the 
probe  around  the  suction  surface  be  computer  controlled. 

The  profile  coordinates,  used  by  a  computerized  milling 
machine  to  manufacture  the  blade  moulds,  were  entered  into 
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the  HP  9845B  computer  and  referenced  to  an  arbitrary  point. 
With  the  probe  support  aligned  with  the  center  blade’s  chord 
line,  the  wire  was  manually  moved  to  a  known  location  in  the 
test  section.  The  traversing  system's  X  and  Y  position 
encoders  were  then  adjusted  to  reflect  this  position  rela¬ 
tive  to  the  same  arbitrary  point.  By  knowing  the  location 
of  the  probe  and  the  blade  surface,  the  computer  could  now 
control  the  movement  of  the  former  in  relation  to  the  lat¬ 


ter  . 


Accuracy  of  the  traverse  was  very  much  dependent  on  the 
accuracy  of  the  Initial  setting  of  the  probe  at  the  known 


location  in  the  test  section.  This  was  estimated  to  be  +/- 
0.005  in  in  both  the  X  and  Y  directions.  The  orientation  of 


the  probe  in  relation  to  the  cascade  during  the  traverses  is 
depicted  in  Fig.  10. 

A  traverse  at  a  particular  chord  location  would  begin 
with  the  computer  positioning  the  probe  wire  at  a  point  0.03 
in  away  from  the  surface  of  the  blade.  Some  vibrations  were 
present  when  the  probe  was  moved  from  one  data  point  to  the 
next.  Though  these  vibrations  dampened  out  quickly,  moving 
the  probe  no  closer  than  0.03  in  from  the  blade  reduced  the 
chances  of  the  wire  contacting  the  surface. 

A  secondary  reason  for  the  offset  starting  position  was 
to  avoid  any  possible  influence  of  the  cooler  blade  surface 
on  the  heat  transfer  characteristics  of  the  probe.  Oka  and 
Kostic  (19:31,33)  found  that  when  a  probe  was  in  close 
proximity  to  a  colder  wall,  the  additional  cooling  effects 
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gave  rise  to  higher  apparent  velocities.  In  their  investi¬ 
gation,  velocities  measured  In  this  region  were  on  the  order 
of  five  to  six  times  greater  than  the  true  value.  The  maxi¬ 
mum  distance  from  the  wall  at  which  this  influence  was  de¬ 
tected  was  0.04  In  for  the  condition  of  zero  velocity  flow. 
In  addition,  the  region  of  Influence  decreased  with  In¬ 
creased  Reynolds  number.  Consequently,  an  offset  distance  of 


0.03  In, 

for 

the 

Reynolds  number 

range  of  this 

study,  was 

felt  to 

be 

well 

clear  of  any 

Influence 

of 

the  wall. 

From 

Its 

starting  position,  the  computer 

would  automa- 

tlcally  move  the  probe  away  from  the  blade  In  increments  of 
0.005  In  for  a  total  of  60  data  points  per  traverse.  At  the 
|  same  time,  the  axial  position  would  be  adjusted  so  as  to 

maintain  a  traverse  normal  to  the  blade  surface.  The  tra- 
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verses  were  depicted  in  Fig.  9.  Once  a  particular  traverse 
was  completed,  the  probe  would  automatically  be  moved  to  the 
next  chordwise  location  and  the  sequence  repeated. 

Unlike  the  wake  survey,  it  was  found  that  consecutive 
readings  of  the  DC  and,  particularly,  the  AC  voltage  output 
varied  considerably  at  each  data  point  close  to  the  blade. 
This  was  thought  to  be  due,  at  least  in  part,  to  the  higher 
frequency  response  of  a  hot  wire  compared  to  a  hot  film 
probe,  requiring  that  some  sort  of  mean  reading  be  used. 
Therefore,  at  each  data  point  in  the  boundary  layer  survey, 
ten  readings  of  the  AC  and  DC  voltage  outputs  were  taken  and 
averaged.  From  a  number  of  trial  runs,  it  was  found  that 
the  arithmetic  mean  of  10  readings  did  not  differ  consider¬ 
ably  from  that  of  30  readings  but  that  the  data-taking  time 
was  significantly  decreased. 

Edge  Velocities .  Outside  of  the  actual  boundary 
layer,  the  measured  velocity  profile  decreased  with  distance 
from  the  blade  due  to  the  presence  of  a  blade  to  blade  nor¬ 
mal  pressure  gradient.  This  complicated  the  determination 
of  the  edge  velocity  and  hence  the  boundary  layer  thickness. 

As  discussed  previously  in  Section  II,  the  measured  vel¬ 
ocity  profile  can  be  assumed  to  be  a  composite  profile  made 
up  of  three  regions.  This  was  mathematically  represented  by 
Eq  (1)  repeated  below: 

Una«a  =  U«l  +  Ukw  -  U«  (  1  ) 

Looking  specifically  at  the  case  right  at  the  wall,  the  no- 
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slip  condition  dictates  that  both  the  measured  velocity  and 
boundary  layer  velocity  equal  zero.  From  the  above  equa¬ 
tion,  UE  must  then  be  equal  to  uiNV  at  the  wall.  There¬ 
fore,  in  order  to  determine  the  edge  velocity  and  subse¬ 
quently  the  boundary  layer  velocity  profile,  one  must  first 
find  a  way  to  extrapolate  the  inviscid  velocity  profile  to 
the  surface  of  the  blade.  Deutsch  and  Zierke  point  out  that 
a  "...  rigorous  approach  does  not  appear  possible  so  one 
must  settle  for  a  self  consistent  attack  which  produces 
plausible  results"  (7:8). 

Ball,  Reid,  and  Schmidt  (2:11-13),  in  their  investiga¬ 
tion  into  the  end-wall  boundary  layer  in  a  two-stage  fan, 
addressed  the  problem  by  using  an  iterative  approach.  They 
began  by  curve  fitting  the  outer  inviscid  velocity  profile 
to  an  assumed  value  for  Ug  .  The  boundary  layer  velocity 
profile  could  then  be  calculated  U3ing  Eq  (1).  Eoundary 
layer  integrations  were  carried  out  to  obtain  the  momentum 
and  displacement  thicknesses  which  in  turn  were  used  to  cal¬ 
culate  the  skin  friction  coefficient  based  on  Ludwieg  and 
Tillman's  relationship  (2:11).  A  subsequent  iteration  of 
Cole's  law  of  the  wake  (2:11)  resulted  in  a  value  for  the 
shear  velocity,  .  A  new  value  for  the  edge  velocity 
could  then  be  obtained  from  the  relationship 

U*  ~  (Cc/2)1'1 

This  value  was  compared  to  the  value  assumed  for  Ug  at  the 
start  of  the  loop.  If  it  was  different,  the  process  was  re- 
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peated  with  the  new  value  until  such  time  the  new  and  old  U 
's  were  approximately  equal. 

This  approach  was  not  possible  in  this  investigation  due 
to  the  unknown  nature  of  the  flow  field  at  distances  closer 
than  0.03  in  from  the  blade  surface.  The  lack  of  data  in 
this  region  precluded  the  determination  of  the  momentum  and 
displacement  thicknesses. 

The  scheme  that  was  employed  in  this  study  was  one  adap¬ 
ted  from  Deutsch  and  Zierke  (7:8-9).  They  found  that  a 
range  of  points  existed  in  the  inviscid  region  through  which 
a  quadratic  curve  fit  extrapolated  to  the  wall  would  yield 
an  edge  velocity  that  was  essentially  invariant.  This  range 
was  represented  by  ninV^2  +/~  where  N IHV  is  the 
number  of  data  points  from  the  maximum  velocity  position  to 
the  measured  point  furthest  from  the  blade.  Williams 
(30:32),  in  his  measurements  of  the  boundary  layer  over  the 
aft  region  of  the  blade,  curve  fitted  ^mv^2  the  outer¬ 
most  points  of  the  measured  velocity  profile.  He  carried 
out  both  a  first  and  second  order  least  squares  fit  and 
chose  the  one  with  the  least  variance. 

Initially,  the  problem  was  approached  in  the  same  manner 
as  Williams.  Difficulties  arose,  however,  when  attempting 
to  duplicate  results  from  one  run  to  the  next.  This  was 
particularly  true  for  traverses  located  near  the  leading 
edge.  In  this  region,  the  inviscid  profile  exhibits  some 
curvature  as  compared  to  the  linear  decay  of  velocity  for 
traverses  further  aft  on  the  blade  (7:6).  The  procedure  was 


therefore  modified  as  follows.  First,  both  linear  and  qua¬ 
dratic  least  squares  fits  were  carried  out  for  three  ranges 


of  outermost  points:  ^IHV^  4,  Ninv/2,  and  3NIMy/4 .  All  of 

these  were  extrapolated  to  the  wall  to  yield  a  maximum  of 
six  values  for  edge  velocity.  Some  second  order  fits  gave 
complex  conjugate  roots  when  extrapolated  and  were  ignored. 
The  remaining  values  for  U^were  reviewed,  and  any  that  were 
significantly  different  from  the  rest  were  deleted.  What 
was  left  was  averaged  to  give  a  mean  value  for  Ug  for  that 
particular  traverse  location. 

Boundary  Laver  Thickness .  Using  Eq  (1),  the  boundary 
layer  velocity  profile  was  calculated  for  each  value  of  edge 
velocity  remaining  from  the  curve  fits.  In  each  case,  the 
boundary  layer  thickness  <5^^  was  found  by  determining  the 
distance  from  the  blade  where 


U*.  -  0.99  Ut 


(20) 
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Note  that  the  particular  edge  velocity  corresponding  to  that 
curve  fit  was  used  and  not  the  mean  value.  The  separate  va¬ 
lues  for  dgL  were  then  averaged  to  give  a  mean  boundary 
layer  thickness  for  that  chordwise  location. 

Loca 1  Turbulence  Intensity.  In  a  manner  similar  to  the 
wake  survey,  the  local  turbulence  intensity  in  the  boundary 


Sr  ^ 

•V  > 1 


layer,  Tug^was  calculated  using 


TUm.  3  Um/Uum 


where  U, 


is  the  Local  time  varying  velocity  measured  3 ■ 


the  local 


the  true  RMS  output  of  the  anemometer  and  U^g^  the  local 
streamwlse  measured  velocity  determined  from  the  DC  output. 

Deutsch  and  Zierke  (7:9)  point  out  that  turbulence  in¬ 
tensities  in  a  classical  equilibrium  boundary  layer  will 
show  a  higher  than  freestream  value  for  distances  y  <,  1.25(5, 
This  relationship  was  used  to  obtain  a  rough  value  for  the 
boundary  layer  thickness  to  compare  with  the  hot  wire  re¬ 
sults. 
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V.  Results  and  Discussion 


Suction  Surface  Pressure  Distribution 

The  static  pressure  measured  at  each  of  the  38  pressure 
tap  locations  on  the  suction  surface  was  reduced  to  pressure 
coefficients  using  Eq  (2).  For  each  blade  configuration, 
the  pressure  coefficients  obtained  from  three  runs  were 
averaged  to  give  a  mean  pressure  profile.  These  profiles 
are  shown  in  Figures  11  through  13.  Though  these  plots  are 
based  on  an  average  of  the  results.  It  was  found  that  the 
profiles  for  a  particular  configuration  were  very  similar 
from  run  to  run.  This  is  an  indication  of  the  consistency  of 
the  flow  through  the  CTF . 

The  profile  for  configuration  1,  Fig.  11,  is  very  sim¬ 
ilar  to  that  obtained  by  Williams  (30:27)  for  his  smooth 
baseline  blade.  The  additional  pressure  taps  used  in  this 
investigation  provide  for  greater  definition  of  the  profile 
particularly  in  the  forward  region  of  the  blade  (seven  taps 
versus  two  used  by  Williams).  The  profile  between  the  25 
and  50  percent  chord  locations  is  relatively  flat  suggesting 
that  the  flow  has  separated.  It  has  not,  however,  com¬ 
pletely  separated  since  the  pressure  coefficient  begins  to 
increase  again  aft  of  the  65  percent  chord  point.  Rather, 
the  flat  region  marks  the  presence  of  a  laminar  separation 
bubble.  The  exact  locations  of  the  points  of  separation, 
transition,  and  turbulent  reattachment  are  difficult  to 


pinpoint  due  to  the  lack  of  the  associated  sharp  discontin¬ 


uities.  It  does  appear  though  that  the  flow  reattaches  by 


about  the  50  percent  chord  point. 


Analyzing  the  pressure  profiles  for  the  two  roughened 


blades.  Figs.  12  and  13,  a  definite  fluctuation  exists  in 


the  pressure  coefficient  between  the  25  and  42  percent  chord 


locations.  This  was  not  apparent  In  Williams  (30:28-30) 


results.  Though  not  definite.  It  is  felt  that  the  cause  of 


this  anomaly  lies  in  the  method  of  applying  the  roughness 


By  Incorporating  the  roughness  in  the  actual  casting  pro¬ 


cess,  It  was  anticipated  that  the  profile  of  the  blade  could 


be  maintained  without  adding  a  joint  or  other  discontinuity 


to  mark  the  edge  of  the  emery  paper.  This,  unfortunately. 


was  not  the  case  as  a  distinct  though  small  ridge  was  appar¬ 


ent  at  about  the  25  percent  chord  point.  This  is  where  the 


roughness  ends,  and  the  regular  surface  of  the  blade  contin¬ 


ues.  The  fluctuations  do  dampen  out,  however,  and  the  pro¬ 


files  aft  of  the  45  percent  chord  location  follow  the  same 


trend  as  for  configuration  1. 


Another  distinct  change  in  the  pressure  profile  with  the 


addition  of  roughness  occurred  In  the  forward  part  of  the 


blade.  The  pressure  coefficient  at  the  forwardmost  tap. 


3.125  percent  chord,  is  significantly  less  than  the  cor¬ 


responding  baseline  value.  For  configuration  2,  there  is 


approximately  a  58  percent  decrease,  while  configuration  3 


experiences  about  a  34  percent  decrease.  The  magnitude  of 


these  jumps  may  have  been  Influenced,  to  a  certain  extent. 


po 


l H 


by  the  location  o£  the  forwardmost  tap  which  coincides  with 
the  transition  from  smooth  to  roughened  surface.  In  both 
cases,  over  a  short  distance  downstream,  the  pressure  in¬ 
creases  to  values  greater  than  that  of  configuration  1. 
This,  of  course,  means  that  the  adverse  pressure  gradient  is 
much  more  severe  for  the  roughened  blades  than  for  the 
smooth  one.  This  can  be  seen  by  the  steep  slope  of  the 
pressure  profile  in  the  forward  region  of  the  blade. 

Due  to  the  fluctuations  present,  it  is  difficult  to 
identify  a  distinct  constant  pressure  plateau  for  configura¬ 
tions  2  and  3  although  one  is  possible  for  the  latter  from 
about  the  38  to  44  percent  chord  points.  Nevertheless,  it 
appears  unlikely  that  the  initial  laminar  boundary  layer, 
with  its  limited  momentum  influx,  could  advance  that  far 
against  such  an  adverse  pressure  gradient  before  separating 
(21:154).  A  more  plausible  interpretation  of  the  results  is 
that  laminar  separation  occurs  shortly  aft  of  the  pressure 
minimum,  where  both  profiles  exhibit  a  local  pressure  maxi¬ 
mum.  These  points  are  at  the  7.8  and  9.4  percent  chord 
locations,  respectively,  for  configurations  2  and  3.  Due  to 
the  finite  spacing  between  pressure  taps,  a  short  constant 
pressure  plateau  may  exist  here.  At  the  very  least,  the 
distinct  variations  In  the  profiles  about  these  points  sug¬ 
gest  that  the  boundary  layer  is  undergoing  some  sort  of 
change,  if  not  a  laminar  separation  bubble,  then  possibly 
just  a  natural  transition  from  a  laminar  to  a  turbulent 
state  without  separation.  For  either  to  occur,  however,  re- 
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quires  the  transition  point  to  move  forward  as  compared  to 
configuration  1.  The  likelihood  of  this  occurring  can  be 
evaluated  by  examining  the  influence  of  surface  roughness  on 
transition. 

When  a  laminar  boundary  layer  becomes  unstable  with  re¬ 
spect  to  small  velocity  disturbances,  it  will  begin  the 
transition  to  the  more  stable  turbulent  mode  (24:334).  The 
distance  over  which  this  transition  takes  place  le  small  and 
often  thought  of  as  a  point,  specifically,  the  transition 
point  (21:159).  The  disturbances  that  affect  the  boundary 
layer  in  this  manner  Include  those  produced  by  surface 
irregularities  such  as  roughness.  If  the  disturbances  are 
large  enough,  the  laminar  boundary  layer  will  become  unsta¬ 
ble  earlier,  and  the  transition  point  will  move  forward.  The 
larger  the  roughness  elements,  the  farther  forward  the  tran¬ 
sition  point  moves  until,  eventually,  its  location  coincides 
with  the  roughness  Itself  (23:536). 

The  above  suggests  that  a  limit  exists  for  surface 
roughness  below  which  no  effect  is  felt  by  the  transition 
point.  This  limit  Is  expressed  as  (24:335) 

^jjT*  <  120  (22) 

where  i/  Is  the  kinematic  viscosity  of  the  flow.  For  a 
typical  inlet  velocity  of  480  ft/sec  and  a  total  temperature 
of  113  0  F,  Eq  (22)  gives  a  critical  roughness  height  of 
13.61  M  m.  Comparing  this  value  to  those  In  Table  I,  both 
configurations  2  and  3  have  equivalent  sand  roughness  values 


that  are  substantially  greater.  Therefore,  it  Is  expected 
that  transition  will  occur  sooner  for  those  two  blade  con¬ 


figurations  as  compared  to  the  smooth  blade  whose  kg  value 
is  less  than  the  critical  height. 

Except  for  a  slight  decrease  with  roughness,  the  pro¬ 
files  aft  of  the  45  percent  chord  for  configurations  2  and  3 
are  essentially  the  same  as  the  baseline  case.  This  rein¬ 
forces  the  idea  that  the  flow  in  this  region  is  also  turbu¬ 
lent.  It  is  possible  for  the  flow  to  be  completely  separa¬ 
ted  aft  of  the  84.375  percent  chord,  but  there  were  no  pres¬ 
sure  taps  located  there  to  verify  that  condition.  It  is  un¬ 
fortunate  that,  due  to  the  blade  mounting  pins,  two  of  these 
gaps  exist  in  the  pressure  data.  This  is  particularly  true 
for  the  forward  one  which  is  close  to  the  leading  edge  of 
the  blade  and  could  possibly  mask  significant  information. 

Boundary  Laver 

Three  series  of  boundary  layer  traverses  were  carried 
out  for  each  configuration  over  a  period  of  two  to  three 
days.  Representative  plots  from  one  of  those  runs  are 
attached  as  Appendices  D,  E,  and  F  for  configurations  1,  2, 
and  3  respectively.  Each  figure  consists  of:  (1)  the  mea¬ 
sured  velocity  profile  showing  the  influence  of  the  normal 
pressure  gradient;  (2)  the  better  of  a  first  or  second  order 
least  squares  fit  of  the  outermost  points;  and 
(3)  the  resultant  boundary  layer  velocity  profile  as  calcu¬ 
lated  from  Eq  (1).  The  inviscid  curve  fit  is  shown  extrapo- 


lated  to  the  wall,  giving  the  value  for  the  edge  velocity 
for  the  particular  run  presented.  For  comparison,  the  mean 
value  for  the  edge  velocity,  based  on  an  average  of  as  many 
as  six  different  curve  fits  per  location,  over  three  sepa¬ 
rate  runs,  is  marked  as  an  E  on  the  velocity  axis. 

Edge  Velocity.  The  mean  edge  velocities  for  each  con¬ 
figuration  are  tabulated  in  Table  II.  Included,  for  compar¬ 
ison,  are  the  corresponding  velocities,  ,  determined  from 
the  pressure  profiles  using  Eq  (19).  The  other  parameters 
listed  will  be  discussed  later. 

The  relatively  good  agreement  between  the  values  of  Ug 
and  indicates  that  the  method  of  obtaining  edge  veloci¬ 

ties  from  the  hot  wire  data  is  a  reasonable  one.  The  best 
agreement  was  obtained  for  those  regions  with  a  flatter 
pressure  profile.  An  example  is  configuration  1  between  the 
25  and  50  percent  chord  locations.  The  largest  differences 
were  at  the  forwardmost  points  which  coincide  with  the  lar¬ 
ger  pressure  gradients.  Data  obtained  by  Deutsch  and  Zierke 
(7:16)  and  Williams  (30:33)  behave  in  a  similar  fashion. 

The  variation  of  the  edge  velocity  along  the  suction 
surface  can  be  seen  in  Fig.  14.  For  this  plot,  the  edge 
velocity  is  non-d imens ional ized  with  respect  to  the  inlet 
velocity.  With  the  exception  of  the  forward  region,  there 
is  little  difference  between  the  three  configurations,  and 
particularly  between  the  two  roughened  ones.  This  is  under¬ 
standable  given  the  closeness  of  the  pressure  profiles.  As 
expected,  the  edge  velocity  decreases  with  distance  along 
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Boundary  Layer  Parameters 


Conf . 


%  Chord 

U«  (ft/uc) 

Ui  (  ft/fK  ) 

6*.  (in) 

Tu  6m.  ( 

4.688 

603.30 

656.61 

0.032 

0.032 

9.375 

592.44 

614.93 

0.038 

0.039 

25.000 

575.94 

577.23 

0.046 

0.048 

29.688 

574.17 

573.27 

0.046 

0.052 

34.375 

570.83 

567.46 

0.046 

0.052 

40.625 

568.37 

564.23 

0.048 

0.048 

45.313 

566.24 

560.45 

0.047 

0.043 

50.000 

562.11 

555.32 

0.048 

0.044 

65.625 

540.01 

520.64 

0.054 

0.055 

70.313 

533.13 

508.93 

0.058 

0.056 

75.000 

528.52 

498.07 

0.066 

0.060 

79.688 

519.72 

487.31 

0.070 

0.062 

84.375 

508.80 

476.27 

0.076 

0.068 

4.688 

580.48 

637.16 

— 

— 

9.375 

573.47 

603.33 

— 

0.031 

25.000 

564.16 

585.68 

0.033 

0.036 

29.688 

559.89 

572.43 

0.033 

0.031 

34.375 

558.59 

572.16 

0.034 

0.039 

40.625 

553.49 

561.47 

0.040 

0.048 

45.313 

551.67 

559.04 

0.043 

0.051 

50.000 

547.64 

550.93 

0.046 

0.056 

65.625 

525.08 

517.38 

0.055 

0.058 

70.313 

517.81 

505.74 

0.059 

0.064 

75.000 

512.05 

495.77 

0.068 

0.060 

79.688 

503.14 

484.57 

0.073 

0.067 

84.375 

495.17 

474.92 

0.079 

0.072 

4.688 

589.81 

667.46 

— 

— 

9.375 

567.71 

574.19 

0.030 

25.000 

560.16 

564.70 

0.037 

0.036 

29.688 

560.43 

575.09 

0.039 

0.045 

34.375 

557.72 

571.67 

0.044 

0.047 

40.625 

553.07 

561.94 

0.050 

0.059 

45.313 

550.19 

557.25 

0.051 

0.063 

50.000 

546.42 

552.69 

0.051 

0.064 

65.625 

523.61 

517.35 

0.062 

0.067 

70.313 

516.92 

506.11 

0.070 

0.068 

75.000 

508.55 

494.46 

0.075 

0.073 

79.688 

500.22 

483.30 

0.080 

0.075 

84.375 

490.32 

475.67 

0.086 

0.080 

the  blade.  It  also  appears  that  edge  velocity  Increases 
with  the  addition  of  roughness  since  the  plot  shows  configu¬ 
ration  3  with  slightly  higher  velocities  than  conf igurat  ion 
2.  This  is  in  keeping  with  the  pressure  profiles  which  show 
configuration  3  having  slightly  less  pressure  over  mo3t  of 
the  surface.  The  increase  in  velocity  is  due  to  the  thick¬ 
ening  of  the  boundary  layer  with  roughness  which  constricts 
the  lnviscid  region  between  the  blades.  In  order  to  main¬ 
tain  the  mass  flow  rate,  the  velocities  in  that  region  must 
increase.  Note,  however,  that  edge  velocities  for  both 
rough  blades  are  lower  than  that  for  the  smooth  blade.  This 
contradiction  is  more  than  the  likely  a  result  of  the  uncer¬ 
tainty  in  the  traversing  mechanism. 

Boundary  Layer  Thickness .  The  boundary  layer  thickness, 
<5bl  /  was  calculated  using  Eq  (20).  Fig.  15  shows  the  vari 
ation  in  boundary  layer  growth  with  roughness.  The  results 
are  also  listed  in  Table  II.  The  boundary  layer  for  config 
uration  1  appears  to  confirm  the  previous  analysis  that 
flow  is  turbulent  over  the  back  half  of  the  blade.  Thi3  can 
be  seen  by  the  distinct  thickening  of  the  boundary  layer  in 
this  region.  It  is  interesting  to  note  that  the  '-onstin' 
pressure  plateau,  associated  with  a  laminar  separation  t  it 
ble,  in  the  pressure  profile  corresponds  to  a  relatively 
flat  boundary  layer  in  Fig.  15. 

The  lack  of  data  points  in  the  forward  reu.  r  :  *  *  *■ 

blade  for  configurations  2  and  3  indicates  1  oca  1  t  m  iar  . 
layer  thicknesses  that  are  less  than  0.03  in,  *  n*»  mi-  cc 
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1 '.stance  the  pr  >be  13  brouqht  from  the  surface.  In  view  of 

trie  Inherent  errors  Involved  In  positioning  the  probe,  how¬ 
ever,  it  is  entirely  possible  for  the  boundary  layer  to  be 
tni  -*er  than  0.03  in  in  this  area.  In  comparison  with  con- 
figurati:*n  1,  the  boundary  layer  growth  shown  for  the  two 
rough  configurations  appears  to  be  that  of  a  developing  tur¬ 
bulent  boundary  layer.  This  would  be  in  keeping  with  the 
1 nte r pretat i  ?n  of  the  corresponding  pressure  profiles  dls- 
.russed  earlier.  As  expected,  configuration  3,  with  Its 

greater  roughness,  has  a  thicker  boundary  layer  than  config- 

ir  it  Ion  Z  However,  measurements  indicate  that  both  rough 

blades  have  -n.nner  boundary  layers  than  configuration  1 
ver  -he  firs-  portion  of  the  blade.  This  region  extends 

*  -  at  u-  'he  37  and  *3  8  percent  chord  points  for  config 
ira-i  ns  i  and  Z  respectively  The  reason  for  this  anomaly 
.*  ir  -.eir  bu-  -.way  be  a  result  ?f  one  or  more  of  the  follow 

'i  1  . - : *  i a  I  probe  positioning  accuracy  estimated  to  be 

•  '  •  ne  s  1-  boh.  the  X  and  Y  directions;  (21  an  ap 

'  -  Tease  i  -  boundary  layer  t  h  i  -  k  r.e  3  s  for  con  f  i  gur  a 

‘  •  i  « e  •  *  -  e  .  ami  -  a  r  separat  i  ;r.  bubble;  and  3  an  a  p 

«„  r  „  ,  ..  -  a  r  ,  q  r  blade's  *  h  :  ■  ■  k  n e  s  -  J  jp  *  lift; 

••  *  ■  *  '4  t  ^  -  r'  *  w r  r  **  *  ^  »•  /  f*  '  .  f 
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slty  in  the  boundary  layer,  as  calculated  by  Eq  (21),  is 
shown  plotted  as  a  function  of  distance  from  the  surface  in 
Appendices  G,  H,  and  I,  corresponding  to  configurations  1, 
2,  and  3,  respectively.  As  mentioned  in  section  IV,  the 
turbulence  intensity  will  increase  from  its  free-stream 
value  as  the  distance  to  the  blade  surface  decreases.  This 
increase  is  first  felt  at  a  distance  of  y  =  1.25  <3 __ .  The 
boundary  layer  thicknesses  obtained  from  this  relationship 
are  tabulated  for  each  configuration  in  Table  II.  The 
agreement  between  thickness  predicted  in  this  manner  and 
that  determined  from  the  velocity  profile  is  inconsistent. 
In  some  cases,  the  two  agree  exactly  while  in  others,  the 
difference  can  be  as  high  as  30  percent.  The  range  is  not 
surprising  when  one  considers  that  the  above  relationship 
for  the  increase  in  turbulence  intensity  assumes  a  classical 
equilibrium  boundary  layer,  i.e.,  one  with  no  pressure  gra¬ 
dient  along  the  surface.  Clearly  the  pressure  profiles  in¬ 
dicate  that  the  suction  surface  boundary  layer  cannot  be 
considered  to  be  in  equilibrium  (7:5,9).  In  addition,  the 
measured  turbulence  intensities  continue  to  decrease  some¬ 
what  outside  of  the  boundary  layer.  Therefore,  unlike 
Oeutsch  and  Zierke's  results  (7),  the  lack  of  a  constant 
tr»e  stream  value  for  turbulence  makes  it  difficult  to  pin¬ 
point  the  location  where  turbulence  intensity  first  increa¬ 
ses  near  the  boundary  layer. 
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Wake  Surveys 

As  detailed  previously,  the  flow  exiting  the  cascade  was 
probed  by  a  two-dimensional  hot  film  sensor  at  four  distan¬ 
ces  downstream  of  the  center  blade's  trailing  edge.  Both 
velocity  and  turbulence  intensity  measurements  were  re¬ 
solved  into  their  respective  X  and  Y  components,  and  the  re¬ 
sultant  vectors  plotted.  The  results  are  collected  in  Ap¬ 
pendices  J,  K,  and  L  corresponding  to  configurations  1,  2, 
and  3,  respectively.  Each  figure  contains  both  velocity  and 
turbulence  intensity  vectors:  the  lighter  lines  for  velo¬ 
city  and  the  darker  ones  for  turbulence  intensity.  The  ori¬ 
gin  of  each  vector  is  the  actual  surveyed  position,  and 
their  magnitudes  can  be  converted  to  their  respective  units 
with  the  scale  factors  indicated.  These  are  150  ft/sec/in 
for  velocity  and  5  percent/in  for  turbulence  intensity. 

The  wake  produced  by  the  boundary  layers  along  the  blade 
is  clearly  discernible.  By  numerically  integrating  the  .  e  - 
locity  profiles,  the  displacement,  momentum,  and  pseudoen 
ergy  thicknesses  can  be  obtained  as  defined  by  Egs  (4) 
through  (6).  These,  along  with  other  wake  characteristics, 
can  then  be  analyzed,  and  the  effects  of  surface  roughness 
determined . 

Wake  Min lmum  Velocity.  The  plot  of  the  downstream  vari 
ation  in  the  velocity  ratio,  vM[N//vO  »  15  shown  in  F.j  1  c 

for  all  three  roughness  conf  igurations  .  For  comparison,  *  we¬ 
nt  Moe ' s  (17)  configurations  are  also  plotted.  His  lata  ; 
for  the  same  blade  type  and  flow  conditions,  tut  wiM,  n 
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pressure  taps  on  the  surface.  One  Is  his  smooth  baseline 

configuration  with  a  measured  arithmetic  average  roughness 
of  0.09  micrometers  (compared  to  0.45  micrometers  in  this 
study),  and  the  other  has  roughness  applied  giving  an  Ra 
value  of  17.9  micrometers.  The  latter  is  comparable  to  con- 
figurition  3  in  this  study  which  has  an  Ra  value  of  18.3 
micrometers.  As  well,  Moe  had  some  results  for  a  measuring 
plane  closer  to  the  trailing  edge. 

In  addition  to  the  data  points,  three  lines  are  shown  in 
the  figure.  The  solid  line  is  a  plot  of  Eq  (9)  with  Lie- 
bleln  and  Roudebush's  constants  of  0.13  and  0.025  for  a  and 
b  ,  respectively.  As  can  be  seen,  it  does,  in  fact,  give 
a  good  representation  of  the  velocity  recovery  in  the  wake 
for  a  cascade  .  A  definite  trend  for  the  effects  of 
roughness  is  also  apparent .  It  seems  that  increasing  sur¬ 
face  roughness  also  Increases,  to  a  certain  extent,  the  dis 
tance  downstream  required  to  recover  the  velocity.  This  is 
not  to  say  that  the  mixing  process  is  slower.  Rather,  with 
a  thicker  boundary  layer  and  corresponding  greater  velocity 
defect  a*"  the  trailing  edge,  the  mixing  process  behind  a 
rougher  blade  takes  a  little  lonqer  to  equalize  the  v^Ij 
lty.  It  Is  Interesting  to  note,  however,  that,  regardless 
t  ‘he  roughness,  appr  >>x  1  nwj  t  e  1  y  9<'  per  ent  ,  t  the  boundary 
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(14:1)  that  the  wake  Is  reenergized  very  quickly  with  the 
majority  of  the  mixing  losses  occurring  within  1/4  to  1/2 
chord  length  downstream  of  the  trailing  edge. 

The  other  two  lines  on  the  figure  are  the  bounding 
curves  in  the  form  of  Eq  (9)  for  the  data  presented.  Con¬ 
stants  for  the  upper  curve  are  a  =  0.1  and  b  =  0.02,  and 
for  the  lower  curve  a  =  0.185  and  b  =  0.037.  In  addi¬ 
tion,  note  the  excellent  agreement  between  configuration  3 
and  Moe's  blade  of  similar  roughness. 

Wake  Form  Factor  ■  Fig.  17  Is  a  plot  of  the  downstream 
variation  of  the  wake  form  factor.  Again,  the  results  from 
Moe's  study  are  included  for  comparison.  The  results  can  be 
related  to  the  variation  of  the  velocity  ratio  through  Eq 
(10'.  These  are  plotted  for  the  same  three  cases  presented 
in  Fig.  16.  The  form  factor  at  the  trailing  edge  for  the 
solid  line  was  determined  by  Llebleln  and  Roudebush  (14:8) 
to  tie  2.8.  The  corresponding  values  for  the  two  limiting 
urves  were  calculated  from  the  data  presented  in  this  stu¬ 
dy  Since  the  upper  curve  in  Fig.  16  closely  approximates 
Moe's  results,  each  of  those  data  points  were  substituted 
inf"  Eg  11  f  i  solve  for  H^»g.  The  average  of  those  results 
iav»-  t  vu  ;  i  e  f  2.19  for  the  trailing  edge  form  factor. 

:  ;  m  i  i  » r  '.  ,  ,  *  he  1  i'.i  points  tor  .-ont  i  jurat  l  on  3  were  used  to 

'  h  t  a  .  <  me  i  n  /  a  '  ■  i  e  of  2.86  for  the  1  o  we  r  -ur  ve  '  r>  trailing 
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The  equation  does  appear  to  predict  the  variation  in  form 
factor  quite  well  for  Moe ' s  smooth  blade.  Its  accuracy  dim¬ 
inishes,  however,  with  roughness  added,  particularly  in  pre¬ 
dicting  values  for  the  first  chord  length  downstream. 

Lieblein  and  Roudebush  (14:9)  observed  that  for  chord 
length  distances  greater  than  0.4,  the  values  of  form  factor 
are  less  than  1.2,  regardless  of  the  value  at  the  trailing 
edge.  This  does  not  appear  to  hold,  however,  beyond  a  cer¬ 
tain  roughness  level  since  the  form  factor  for  configuration 
3  is  greater  than  1.2  at  a  chord  length  distance  of  0.6.  In 
all  cases  though,  the  wake  form  factor  decreases  rapidly 
just  aft  of  the  trailing  edge  and  appears  to  approach  a  val¬ 
ue  of  1  asymptotically  far  downstream  as  predicted  by  Eq 
(10)  . 

Wake  Full  Thickness .  The  full  thickness  ratio,  6/c  ,  is 
shown  plotted  against  downstream  distance  in  Fig.  18.  As 
expected,  blades  with  greater  roughness  have  thicker  wakes. 
The  mixing  process  that  occurs  in  the  wake  is  shown  by  the 
Increase  in  the  full  thickness  ratio  with  downstream  dis¬ 
tance.  Once  again,  note  the  excellent  agreement  between  the 
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factor,  data  for  each  configuration  was  used  to  calculate  a 

mean  value  of  (<3/c  Vr  £  or  that  particular  roughness  configu¬ 
ration.  The  line  shown  is  a  plot  of  Eq  (11)  directly. 

The  relation,  based  on  isolated  airfoil  data,  does  in¬ 
deed  apply  to  cascades  of  smooth  blades.  This  is  shown  by 
the  excellent  agreement  with  Moe ' s  smooth  blade  and,  to  cer¬ 
tain  extent,  configuration  l's  data.  The  agreement  gets 
progressively  worse,  however,  with  increased  values  of 
roughness  with  one  exception.  At  a  chord  length  distance  of 
1.125,  the  full  thickness  ratio  for  all  five  configurations 
presented  falls  very  close  to  the  line. 

Wake  Momentum  Th  lckness  .  The  plot  of  *roe  momentum 
thickness  ratio,  d  2  c  >  versus  downstream  distance  is  shown 
in  Fig.  20.  For  the  most  part,  the  moment  in-  ‘hie* ness  .5 
essentially  constant  for  each  configuration.  This  agrees 
with  Liebieln  and  Foudebush  14:1  c  11  who  the  r 0  *  :  •  a  '.  y 
showed  that  for  x  c  ■  0.2  ,  little  change  •••'*.•  jr  ?  .  r.  *  2  ', 
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DISTANCE  DOWNSTREAM  OF  TRAILING 


From  Fig.  20,  it  is  clear  that  surface  roughness  increa¬ 
ses  momentum  thickness  in  the  wake  with  the  largest  jump  oc¬ 
curring  between  Moe's  smooth  blade  and  configuration  1.  This 
same  kind  of  jump  is  detectable  in  all  the  wake  parameters 
examined  previously,  leading  one  to  conclude  that  an  initial 
Increase  in  roughness  produces  the  greatest  change  in  wake 
character istics .  Moe  (17:47)  reported  similar  observations 
In  hi3  study  with  regards  to  the  effects  of  roughness  on  the 
total  pressure  loss  coefficient. 

Total  Pressure  Loss .  Using  the  procedure  outlined  in 
Section  IV,  the  total  pressure  loss  coefficient,  0  ,  was  de¬ 
termined  for  each  downstream  measuring  station  from  Eq  (3). 
Tee  results  are  tabulated  in  Table  III.  Included  in  the 
.'.sting  are  total  pressure  loss  coefficients  calculated  from 
obtained  by  Moe  (17)  and  Williams  (30).  The  latter  was 
♦  !  !e  i  since  that  blade's  roughness  level  was  also  similar  to 
:  .  i  .Mi  ion  3  . 

.  f  . -ant  discrepancies  are  evident  in  the  results 
One  example  is  the  variation  in  B  among  the 
i  1“  in f i gurat  i  ons  of  similar  roughness,  i.e.,  con- 
•  .  ►>  r  3,  Moe's  configuration  number  13,  and  Wil- 

r  i*  .  n  number  3.  Values  of  B  calculated  in 
«  i:-,  n  average,  a ppr ox i ma t e 1 y  92  and  73 

t  ••  *■  f  Moe  and  Williams,  respective- 

;  e  n  t  s  tor  configuration  1 , 

.  , 
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ir^i^r  than  Moe's 


Sunnarv  of  Total  Press' 


Conf  . 


Poulin 


Poulin 


Poulin 


Vi Ilians 

Poulin 

(THEORY) 

Poulin 

(THEORY) 

Poulin 

(THEORY) 


Ra,  M 


0.45 


12.1 


18.3 


0.09 


17.9 


19.5 


0.45 


12.1 


18.3 


.125  0.625  1.125  1.625  2.125 


0.0641 


0.0729 


0.0808 


0.0403  0.0423 


0.0211 


0.0780 


0.0895 


0.0398 


0.0329  0.0455  0.0457 


0.0559  0.0516 


0.0212 


0.0659 


0.0770  0.0750 


0.0877 


0.0409 


0.0461 


0.0484 


0.0213 


0.0663 


0 . 0770 


0.0896 


0.0381 


0.0435 


0.0465 


0.0235  0.0241  0.0246  0.0250 


0.0265  0.0268 


the  power  velocity  assumption  and  Eq  (14)  are  considered  to 
be  valid. 


The  theoretical  total  pressure  loss  coefficients,  calcu¬ 
lated  from  Eq  (12)  for  the  three  blade  configurations  used 
in  this  investigation,  are  listed  in  Table  III.  It  appears 
that  d iscrepanc ies  also  exist  between  these  theoretical  re¬ 
sults  and  actual  values  obtained  experimentally.  The  sig¬ 
nificant  spread  between  the  two  can  be  appreciated  from  Fig. 
22.  Nevertheless,  a  number  of  qualitative  observations  can 
be  made.  Firstly,  loss  in  total  pressure  increases  with 
downstream  distance  from  the  trailing  edge.  This  can  be 
seen  as  an  increase  in  Q  for  the  theoretical  curves  in 
Fig.  22.  The  increase  is  very  slight  indicating  that  the 
majority  of  the  mixing  losses  occurred  at  distances  closer 
to  the  trailing  edge.  This  is  in  keeping  with  the  analysis 
of  the  wake  velocity  recovery  made  earlier.  It  can  also  be 
verified  by  calculating  0  for  the  case  of  complete  mixing 
using  Eq  (15).  The  equation  was  applied  at  the  four  measur¬ 
ing  stations  for  each  configuration  and  the  results  aver¬ 
aged.  It  yielded  values  for  of  0.0215,  0.0246,  and 

0 . 0  2  ’’  4 ,  respectively,  for  configurations  1,  2,  and  3.  These 

ire  not  much  different  than  the  corresponding  values  for  □ 
obtained  a*  each  measuring  station.  Except,  tor  a  sudden  and 
.  n  e  *  p  ,  i  .  r .  e  <  '  .  rrg.  1 1  i  -  h  ■ .  r  d  length  5  ;  s  ’  a  n  ••  f  1  .  1  1  ■  ,  an 


Both  the  theoretical  and  experimental  curves  1  nd  1  ate 
that  an  increase  in  surface  rouqhness  leads  to  an  increase 
in  total  pressure  loss.  This  is  not  surprising  ir.  view  of 
the  thicker  boundary  layer  near  the  trailing  edge  and  qreat 
er  velocity  defects  in  the  wake  associated  with  greater 
roughness  levels.  It  is  Interesting  to  note  that,  except 
for  the  anomaly  at  x/c  =  1.125,  the  percentage  increase  in  0 
1 3  comparable  between  theory  and  experiment.  For  example,  at 
a  chord  length  distance  of  0.625,  there  is  approximately  a 
25  percent  increase  In  total  pressure  loss  between  configu¬ 
rations  1  and  3,  measured  experimentally  or  theoretically. 
Farther  downstream  at  x/c  =  2.125,  the  increase  i3  35  per¬ 
cent  experimental  and  31  percent  theoretical.  Thi3  agree¬ 
ment  indicates  that  even  though  some  questions  exist  with 
the  quantitative  results,  the  method  does  provide  reasonably 
good  qualitative  information. 

Loca 1  Turbulence  I nten3 1 ty .  The  effects  of  surface 
roughness  on  the  local  turbulence  intensity,  Tu  ,  in  the 
wake  can  be  surmised  from  the  plots  in  Appendices  J  through 
L.  It  appears  that  the  initial  application  of  roughness, 
configuration  2,  causes  the  freestream  turbulenc  to  approxi¬ 
mately  double  compared  to  configuration  1.  No  further  in¬ 
crease  in  freestream  turbulence  is  discernable  with  addi¬ 
tional  rouqtnes3,  as  represented  by  configuration  3.  In 
fact,  a  very  slight  decrease  in  Tu  ,  on  the  order  of  one  or 
wo  tenths  of  a  percent,  is  observed  from  configuration  2  to 
’.  other  effects  if  increasing  rouqhness  include  an  in 
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•reds?  In  the  affected  wake  area  and  a  larger  3pread  of  peak 
ti  rctuatlons.  These  general  trends  were  also  ooserved  by 

Mop  ■.  i  7  i  . 

The  most  distinguishing  features  of  the  wake  turbulence 
intensity  profiles  are  the  large  fluctuations  emanating  pri¬ 
marily  from  the  suction  surface  of  the  blade  (lower  cross 
stream  location).  The  pressure  surface,  on  the  other  hand, 
though  still  contributing  to  the  affected  wake  area,  has 
fluctuations  that  are  smaller  in  magnitude  and  essentially 
streamwise  in  orientation.  Referring  to  Fig.  105,  the  wake 
profiles  for  con f 1 gurat ion  1,  It  Is  interesting  to  note  the 
magnitude  and  orientation  of  the  turbulence  intensity  fluc¬ 
tuations  from  the  suction  surface  of  a  blade  that  has  no  ad¬ 
ded  roughness.  The  cause  of  the  peaks  could  be  Interpreted 
to  be  the  38  pressure  taps  drilled  in  the  suction  surface 
except  that  the  same  trend  can  be  round  in  Moe ’ s  profile 
(17:78)  for  his  smooth  blade  with  no  pressure  taps.  It  ap¬ 
pears,  therefore,  that  Fig.  105  is  representative  of  the 
turbulence  intensities  found  over  a  smooth  blade  of  this 
type  under  these  conditions.  The  same  Tu  "footprint"  is 
observed  in  the  profiles  of  configurations  2  and  3  at  the 
same  downstream  distance,  Flg3.  109  and  113,  respectively. 
Superimposed  on  thi3  "footprint"  are  the  effects  attribu¬ 
table  to  the  addition  of  surface  roughness.  For  configura 
tlon  2,  the  Tu  vectors  show  a  larger  spread  in  orientation 
but  no  increase  In  peak  values.  Configuration  3,  on  the 
other  hand,  has  both  a  laroer  spread  arid  greater  magnitudes. 
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It  would  be  of  interest  to  see  if  the  trends  observed  in 
the  wake  turbulence  intensity  profiles  could  be  traced  back 
to  the  turbulence  intensity  profiles  along  the  suction  sur¬ 
face  .  Reference  is  made  to  Figs.  78,  91,  and  104  which  de¬ 
pict  these  profiles  at  the  84.375  percent  chord  location  for 
conf 1 gurat  1  ons  1,  2,  and  3,  respectively.  It  appears  from 
these  figures  that  the  peak  value  of  turbulence  intensity 
does  Indeed  remain  fairly  constant  from  configurations  1  to 
2  before  increasing  with  configuration  3.  In  addition,  the 
variation  of  freestream  turbulence  with  roughness,  observed 
in  the  wake,  can  also  be  seen  on  the  suction  surface.  That 
is,  the  value  for  freestream  turbulence  increases  from  con¬ 
figuration  1  to  2  before  decreasing  slightly  with  configura¬ 
tion  3.  Note  that  these  are  qualitative  observations  only. 
Turbulence  intensities  measured  along  the  suction  surface 
are  significantly  larger  than  what  was  found  in  the  wake. 
One  reason  for  this  may  be  that  a  hot  wire  was  used  for  the 
boundary  layer  measurements  while  a  hot  film  was  used  in  the 
wake.  Another  contributing  factor  could  be  that  the  large 
fluctuations  dampen  out  somewhat  prior  to  the  first  measur¬ 
ing  station  in  the  wake,  located  1.25  in  downstream  of  the 
trailing  edge.  Note  the  large  amount  of  damping  that  occurs 
between  that  station  and  the  next  one  at  2.25  in.  Similar 
damping  would,  therefore,  be  anticipated  immediately  aft  of 
the  trailing  edge.  Thi3  premise  is  reinforced  by  reviewing 
Mop ' s  (17)  plots  at  a  downstream  location  of  0.25  in.  Each 
of  his  configurations  show  significant,  damp  i  ng  of  *  he  peak 


values  o£  Tu  from  that  station  to  the  next  one  at  1.25  In. 

A  final  note  on  the  wake  turbulence  intensity  profiles 
has  to  do  with  the  downstream  variation.  Regardless  of  sur¬ 
face  roughness,  the  large  peaks  are  found  to  be  completely 
damped  out  by  a  chord  length  distance  of  2.125.  Over  the 
same  distance,  however,  the  affected  area  in  the  wake  under¬ 
goes  little  change  as  compared  to  the  velocity  profile.  This 
suggests  that  the  turbulence  generated  by  the  blade  will 
still  be  felt  at  a  distance  downstream  where  the  velocity 
distribution  is  once  again  uniform. 
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VI .  Conclus Ions  and  Recommendations 

Conclus  ions 

The  major  conclusions  of  this  investigation  are: 

1.  The  transition  point  of  the  boundary  layer  along  the 
suction  surface  of  the  blade  moves  forward  with  the 
addition  of  roughness. 

2.  The  boundary  layer  edge  velocity  as  well  as  the 
boundary  layer  thickness  near  the  trailing  edge  in¬ 
crease  with  roughness. 

3.  The  addition  of  roughness  increases  wake  full  thick¬ 
ness,  momentum  thickness,  and  form  factor  while  de¬ 
creasing  the  velocity  recovery  immediately  down¬ 
stream  of  the  trailing  edge. 

4.  The  total  pressure  loss  coefficient  Increases  with 
roughness  as  well  as  distance  downstream. 

5.  Wake  characteristics  are  most  sensitive  to  an  ini¬ 
tial  Increase  in  roughness. 

6.  Roughness  effects  on  wake  turbulence  intensity  in¬ 
clude  an  increase  in  free-stream  turbulence,  affec¬ 
ted  wake  area,  and  peak  fluctuations.  These  trends 
can  be  predicted  from  suction  surface  turbulence 
intensity  profiles  near  the  trailing  edge. 

In  addition  to  the  above  noted  roughness  effects,  the 
following  general  observations  are  made: 

1.  Laminar  separation  with  turbulent  reattachment 
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appears  to  occur  on  the  suction  surface  ut  the 
smooth  blade  . 

2.  The  method  used  to  determine  the  boundary  layer  edge 
velocity  is  a  reasonable  one,  giving  results  com¬ 
parable  to  that  obtained  from  the  pressure  profile. 

3.  A  reasonable  qualitative  assessment  of  the  boundary 
layer  growth  can  be  interpreted  from  the  correspon¬ 
ding  turbulence  intensity  profiles. 

4.  The  relationship  between  pseudoenergy  thickness  and 
wake  form  factor,  developed  by  assuming  a  power 
profile  for  the  wake  velocity  profile,  is  a  satis¬ 
factory  one,  regardless  of  roughness  levels,  for 
form  factors  up  to  1.2. 

5.  Downstream  variation  of  wake  momentum  thickness  ra¬ 
tio  and  increase  in  full  thickness  ratio  for  cas¬ 
cade  flow  have  similar  trends  to  that  of  flow  over 
an  isolated  airfoil. 

6.  Majority  of  the  changes  in  the  wake  characteristics 
and  mixing  losses  occur  within  half  a  chord  length 
distance  downstream  of  the  trailing  edge. 

7.  Turbulence  generated  by  the  blade  will  continue  to 
be  felt  in  the  wake  even  after  the  velocity  distri¬ 
bution  is  once  again  uniform. 

Recommendat i on5 

As  previously  suggested  by  Moe  (17)  and  Williams  (30), 
the  effects  of  varying  inlet  turbulence  on  the  performance 
of  the  cascade  would  be  of  interest.  One  obiertive  should 
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Though  the  boundary  layer  along  the  suction  surface  of 
the  blade  was  investigated,  and  the  procedure  successfully 
automated,  the  information  obtained  was  deficient  in  that  no 
data  were  available  closer  than  0.03  in  from  the  surface. 
The  primary  reason  for  this  was  concern  in  contacting  the 
blade  surface  with  the  probe.  This  would  be  alleviated  by 
switching  from  a  hot  wire  to  a  laser  anemometer  system.  Of 
course,  any  benefits  gained  would  have  to  be  weighed  against 
the  increased  complexity  of  the  system  and  required  modifi¬ 
cations  to  the  test  section. 

Finally,  discrepancies  in  the  total  pressure  loss  coef¬ 
ficients  obtained  in  this  and  the  previous  two  investiga¬ 
tions  suggest  that  a  problem  exists  in  the  data  reduction 
procedure.  The  difficulty  is  further  emphasised  by  the  dis¬ 
parity  between  theoretical  and  experimental  loss  coeffi¬ 
cients  calculated  in  this  study  from  the  same  wake  velocity 
profiles.  Consequently,  the  procedures  and  software  used 
should  be  reviewed  in  an  attempt  to  either  confirm  the  re¬ 
sults  or  identify  the  problem. 
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APPENDIX  A:  Component 

Listing 

R 

Component 

TvDe/Model  # 

/- 

Pressure  Transducers 

r." 

Tank  Total  Pressure 

Statham  PM60TC 

Test  Section  Inlet  Static 

Statham  P6TC 

■ 

Test  Section  Exit  Static 

Statham  P6TC 

A 

Ambient  Pressure 

CEC  4-326 

Bridge  Balance 

CEC  type  8-108 

r% 

DC  Power  Supply 

HP  6205C 

.N 

Scanlvalve  System 

*T 

Pressure  Transducer 

PDCR  2 3D 

Scanivalve 

48S9-3003 

» . 

Controller 

CTLR  2/S2-S6 

.  % 

Scanner  Position  Display 

J102/J104 

v\ 

& 

Thermocouples 

Four  copper-constantan  types 

to  measure  ambient  temperature, 
tank  total  temperature,  and 
calibrator  chamber  top  and 

Omega  T-type 

bottom  total  temperature 

1 

Traversing  Mechanism 

Motors  (2) 

North  American  P-. 

v;; 

Controls  K  ^  2 ^  L 

Encoder  Transducers  (2) 

Astrosystems  MT ~ 

■ 

Hot  Wire/Film  Anemometer  System 

£ 

Anemometers  (3) 

TS I  Mode  1  1  " c- 

Monitor  and  Power  Supply 

TSI  Model  1 

Osc  1 1 loscope 

B&K  Mode  1  1 •  '  a 

•  > 

v**. 

Boundary  Layer  Hot  Wire 

TSI  Mode  .  ’ 

k“ 

X-con£ igurat ion  Hot  Film 

TSI  Mode:  J 

Boundary  Layer  Probe  Support 

TSI  Mod* . 

F: 

X-conf igurat ion  Probe  Support 

TSI  M  *  »■> 

.■* 

Calibrator  (modified) 

”**  *  M  •  , 

Transformer  (to  heat  air 

£ 

for  calibrator' 

V  4  ‘  ♦ 

k. 

Data  Acquisition  System 

•  ’f 

Computer 

*  '. 

Disk  Dr i ves  ( 2  ! 

to 

Channe 1  Scanner 

Digital  ”•>  1  r  me  *  *  r 

Printer 

P 1  otter 

§ 
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APPENDIX  B:  Hot  Wlre/Fllm  Calibration 

The  procedure  used  to  calibrate  the  hot  wires  and  films 
in  this  investigation  was  one  developed  by  Rivir  and  Vonada 
(29)  to  take  into  account  elevated  temperatures.  This  was 
necessary  since  the  flow  exiting  the  CTF  did  so  at  tempera¬ 
tures  typically  30  to  40  °F  above  ambient.  In  addition  the 
temperature  at  which  the  system  stabilized  for  a  particular 
run  would  vary  from  day  to  day  depending  on  weather  condi¬ 
tions.  For  the  remainder  of  this  appendix,  the  term  "wire" 
will  be  used  to  denote  both  the  hot  wire  and  hot  film  sen¬ 
sors  used  in  the  study. 

Hot  wire  anemometry  is  basically  a  heat  transfer  pro¬ 
blem.  Therefore,  of  primary  interest  in  the  calibration  is 
the  variation  of  Nusselt  number  with  Reynolds  number.  King 
was  one  of  the  first  to  study  this  and  developed  the  follow¬ 
ing  semi-empirical  equation  known  as  King's  law: 

Nu  *  A  ♦  BRe"  (23) 

where  Nu  is  the  Nusselt  number,  A  and  B  the  intercept  and 
slope  of  the  calibration  curve  respectively.  Re  the  Reynolds 
number,  and  n  an  exponent  applicable  for  a  particular  range 
of  Reynolds  number.  King  used  a  value  of  0.5  for  n  (5:114). 

Heat  is  transferred  out  of  the  wire  by  radiation,  buo/- 
ant  convection,  conduction  to  the  end  supports,  and  forced 
convection  by  the  fluid  flow.  Of  these,  the  first  two  are 


neglected  since  radiation  heat  loss  amounts  to  only  about 
about  0.1  percent  of  the  electrical  Input  to  the  wire  while 
buoyant  convection  is  significant  only  at  very  low  speeds 
(5:109).  The  effect  of  end  losses  due  to  conduction  to  the 
supports  and  convection  at  zero  velocity  are  contained  in 
the  quantity  A  in  Eq  (23)  (9:177). 

Referring  to  King's  law  Eq  (23),  Bradshaw  points  out 
that  correlations  of  this  type  "...  are  applicable  only  over 
a  limited  range  of  Reynolds  number  and  there  is  no  physical 
reason  why  a  single  exponent  should  do  for  an  indefinitely 
large  range"  (5:114).  For  this  reason,  the  calibration 
scheme  employed  the  following  values  for  n  chosen  by  Collis 
and  Williams  (6:370) 

n  =  0.45  for  the  range  0.02  <  Re  <  44 
n  =  0.51  for  the  range  44  <  Re  <  140 

"The  temperature  difference  between  the  wire  and  the 
fluid  is  usually  large  enough  for  the  choice  of  reference 
temperature  for  the  fluid  properties  to  be  in  doubt" 
(5:115).  Collis  and  Williams  (6:370)  recognized  this  pro¬ 
blem  and  chose  to  evaluate  fluid  properties  at  a  mean  tem¬ 
perature  defined  to  be  TM  =  (T^  +  Tq)/2  where  T y  is  the 
temperature  of  the  wire  and  T0  the  temperature  of  the  free 
stream.  Their  final  form  of  the  King's  law  equation,  which 
includes  a  so-called  "temperature  loading"  factor  applied  to 
the  Nusselt  number,  reads  as  follows: 


Nu(T»/To)" 


A  +  BRe” 


(24) 


with  the  exponent  m  having  a  value  of  -0.17  (6:369).  The  fi¬ 
nal  result  of  the  calibration  procedure  outlined  below  is  a 
curve  of  this  form. 

Calibration  Procedure 

The  first  step  in  calibrating  a  wire  was  to  determine 
how  its  resistance  varied  with  temperature.  This  was  accom¬ 
plished  by  placing  the  sensor  in  a  low  velocity  flow  main¬ 
tained  at  a  particular  temperature  and  measuring  its  cold 
resistance.  This  was  repeated  for  different  temperature 
settings  in  the  range  of  interest  which  for  this  study  was 
100  to  125  o  F.  By  extrapolating  this  data,  the  required 
operating  resistance  needed  to  yield  a  wire  temperature  of 
500  OF  was  calculated. 

An  important  parameter  in  the  calculation  of  heat  trans¬ 
fer  is  the  adiabatic  wall  temperature  T^w,  also  known  as  the 
recovery  temperature.  It  is  the  equilibrium  temperature  at¬ 
tained  by  a  surface  exposed  to  a  flow  when  there  is  no  heat 
transfer  at  the  surface.  The  adiabatic  wall  temperature  for 
the  wire  was  determined  by  setting  it  in  the  calibrator  flow 
at  the  expected  test  velocity  and  temperature  and  measuring 
its  cold  resistance.  The  expected  test  velocities,  based  on 
previous  investigations,  were  450  ft/sec  for  the  wake  survey 
and  500  ft/sec  for  the  boundary  layer  survey,  both  at  tem¬ 
peratures  of  113  °F.  Since  the  wire  was  not  being  heated  by 
the  anemometer,  it  was  assumed  that  it  would  quickly  reach  a 
temperature  at  equilibrium  with  the  flow,  and  the  condition 


o£  no  heat  transfer  at  the  wall  would  be  attained.  The 
equation  for  the  adiabatic  wall  temperature  is  given  by 

T*.  -  T.  +  <V*V2c,)re  (25) 

where  rc  is  a  nondimens ional  fluid  temperature  known  as 
the  recovery  factor  (11:297). 

Applying  the  isentropic  relationship  between  total  and 
static  temperatures  yields  the  following  equation  for  the 
recovery  factor: 

r«  -  (T«,-T,)/<Tt.-T.)  (26) 

where  TTq  is  the  total  temperature  of  the  flow.  This  val¬ 
ue  for  the  recovery  factor  was  assumed  to  be  constant  and 
was  used  repeatedly  in  the  calibration  as  well  as  data  re¬ 
duction  to  calculate  the  adiabatic  wall  temperature. 

The  question  is  again  raised  as  to  what  temperature 
should  the  fluid  properties  be  evaluated  at.  Due  to  the 
relatively  high  velocities  to  be  measured,  the  reference 
temperature  was  calculated  as  follows: 

Th  -  (T*-T.)/2  ♦  0 . 22(T«m-T«)  (27) 

This  equation,  formulated  by  Eckert  assuming  constant  speci¬ 
fic  heat,  yields  solutions  for  high  velocity,  laminar  boun¬ 
dary  layers  that  are  applicable  over  a  wide  range  of  free- 
stream  and  surface  temperatures  (11:304).  Note  that  in  Col- 
lis  and  Williams  investigation  (6:361),  the  maximum  velocity 
used  was  140  ft/sec.  Applying  Eq .  (25)  at  this  speed  yields 


a  difference  of  only  1.63  °F  between  the  adiabatic  wall  and 
free-stream  flow  temperatures,  assuming  a  Prandtl  number  of 
unity.  Consequently,  for  the  range  of  velocities  studied  by 
Collis  and  Williams,  Tj ^  *  Tq  and  Eq  (27)  reduces  to  the 
relationship  that  they  employed  for  calculating  the  refer¬ 
ence  mean  temperature. 

The  data  points  for  the  actual  calibration  curves  were 
obtained  by  measuring  the  output  voltage  of  the  anemometer 
for  15  to  20  velocities  in  the  range  expected.  These  were 
measured  at  three  fluid  temperatures:  100,  113,  and  125  °F . 
In  addition  to  the  voltage,  the  acquisition  system  sampled 
the  calibrator  chamber  total  pressure,  inlet  and  exit  total 
temperatures,  and  the  ambient  pressure  at  each  data  point. 
With  these  values,  the  velocity  of  the  flow  was  determined 
along  with  the  adiabatic  wall  temperature  using  Eq  (25). 

The  Reynolds  number  plotted  on  the  calibration  curves 
was  based  on  the  diameter  of  the  wire: 

Re  -  <*V«d)/M  (28) 

where  Vc  is  the  calibrator  velocity,  d  the  diameter  of 

the  wire,  and  p  the  dynamic  viscosity.  The  fluid  proper¬ 
ties  were  evaluated  at  the  reference  temperature  calculated 
from  Eq  (27).  The  Reynolds  number  based  on  the  data  taken  at 
the  expected  test  velocity  determined  the  value  of  the  expo¬ 
nent  n  from  Collis  and  Williams  criteria. 

The  final  parameter  to  be  determined  is  the  Nusselt  num¬ 
ber.  In  a  constant  temperature  anemometer  system,  the  wire 


is  maintained  at  a  specific  temperature  by  a  feedback  con¬ 
trolled  bridge  circuit.  At  equilibrium,  an  energy  balance 
of  the  wire  dictates  that  the  electrical  power  generated  by 
the  anemometer  to  maintain  the  constant  temperature  must 
equal  the  heat  dissipated  by  the  wire  to  the  fluid  over 
time.  This  can  be  expressed  as 

R*I*  »  hA(T*-T.)  (29) 

where  Ry  is  the  resistance  of  the  wire,  I  the  current 
through  the  wire,  h  the  convection  heat  transfer  coeffi¬ 
cient,  and  A  the  surface  area  of  the  wire. 

Substituting  Ey  /Ry  for  the  current,  where  Ey  is  the 
voltage  across  the  wire,  and  rewriting  the  heat  transfer 
area  as  the  product  of  the  wire's  circumference  and  its 
length  L,  yields 

E-VR.  -  h*dL(T«-T.)  (30) 

Dividing  through  by  the  thermal  conductivity  of  the  flow  and 
rearranging 


The  expression  hd/k  is  termed  the  Nusselt  number,  based  on 
the  diameter  of  the  wire. 

From  the  schematic  of  the  bridge  circuit.  Fig.  23,  it 


can  be  shown  that 


voltmeter 


Var  table 

Rec 1  stance 
Decade 


Fig.  23.  Simplified  Schematic  of  the  Anemometer  Bridge  Circuit 


where  3o  is  the  voltage  output  of  the  anemometer  system  and 
ER3  the  voltage  drop  across  the  resistor  R3  .  Replacing  ER. 
with  IR3  and  introducing  I  =  Eq  /(Rj  +Ry  )  yields 


Be  *  Bel  Re/  (Re+Re)  ] 


(32) 


Substituting  Eq  (32)  into  Eq  (31)  results  in  the  final  form 
for  the  Nusselt  number  as  a  function  of  the  anemometer  out¬ 
put  voltage: 


B«*Re 


ITU  \  K|TR«  f  - 1  l  e-Te ) 


(33) 


The  calibration  curves  themselves  are  plots  of  the  Nus¬ 
selt  number  times  a  temperature  loading  factor  raised  to  an 
exponent  m  versus  Reynolds  number  raised  to  the  exponent  n. 
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An  example  is  depicted  In  Fig.  24.  In  this  case,  the  three 
separate  curves,  corresponding  to  the  three  temperatures  at 
which  the  calibrations  were  carried  out,  are  clearly  discer¬ 
nible  . 

The  final  step  of  the  calibration  procedure  is  for  the 
computer  to  iterate  to  determine  the  value  of  the  exponent  m 
for  which  a  linear  regression  of  all  the  data  points  yields 
the  smallest  variance.  The  slope  and  intercept  of  that 
curve  fit,  which  correspond  to  the  quantities  B  and  A  in 
Eq  (24),  are  saved  on  magnetic  disks  for  subsequent  use  in 
reducing  anemometer  output  voltages,  obtained  In  experimen¬ 
tal  runs,  into  velocities. 
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Calibration  Curves  at  Three  Temperatures 


APPENDIX  C:  Roughness  De £  1 n 1 1 1 ons 

"The  quality  of  the  surface  is  defined  as  the  height  of 
the  peaks,  the  structure  of  the  roughness  and  the  waviness. 
In  workshop  practice,  however,  only  a  roughness  value  is 
prescribed"  (4:30).  The  most  common  of  these  roughness 
parameters  Is  the  arithmetic  average  roughness,  Ra,  defined 
to  be  the  average  value  of  the  deviations  from  the  center- 
line  over  a  particular  sampling  length  (20:5).  Schaffler 
(22:10)  expressed  this  mathematically  as 

fL 

Ra  “  r  jlyl  dx  (34) 

<r 

where  L  is  the  roughness  sample  length,  and  x  and  y  the 
roughness  length  and  height  coordinates,  respectively,  as 
shown  below: 


Fig.  25.  Average  Roughness,  Ra 

The  traditional  way  of  specifying  roughness  in  fluid 
mechanics  is  in  terms  of  an  equivalent  sand  roughness,  kg  . 
By  correlating  the  results  of  two  previous  studies,  Koch  and 
Smith  (12:424)  arrived  at  the  following  equation  to  relate 


* 


equivalent  sand  roughness  to  profllometer  readings: 


f 


£ 


t 

I 


k,  >  6.2  Ra  (35) 

Schaffler  (22:10)  points  out,  however,  that  there  are 
some  problems  in  using  Ra  or  kg  to  quantify  the  tech¬ 
nical  roughness  associated  with  various  manufacturing  pro¬ 
cesses.  For  example,  equivalent  sand  roughness  is  based  on 
tightly  packed  grains  of  equal  size  while  technical  rough¬ 
ness  consists  of  a  large  range  of  peak  heights  and  surface 
waviness.  The  calculation  of  the  arithmetic  average  value 
of  roughness,  on  the  other  hand,  includes  small  roughness 
elements,  the  effects  of  which  are  insignificant  in  defining 
the  hydrodynamic  character istic  of  the  surface.  Consequent¬ 
ly,  Schaffler  proposes  a  different  roughness  parameter,  k  , 
which  describes  the  peaks  rather  than  the  average  value  of 
roughness.  "Roughness  k  is  defined  as  the  difference  be¬ 
tween  the  arithmetic  averages  of  the  ten  highest  peaks  and 
the  ten  deepest  grooves  which  exist  per  millimeter  length" 
(22:10)  or 


a 

-v 


k  ■  y*M«  -  y«Mm 

Figure  26  is  an  example  of  a  roughness  trace  with  the 
roughness  height  k  depicted.  By  measuring  the  surface 
roughness  of  a  large  number  of  blades,  he  obtained  the 
following  simple  correlation  between  k  and  Ra : 

k  =  8 . 9  Ra  ( 36 ) 
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Fig,  33.  Boundary  Layer  Velocity  Profiles,  Conf.#l  at  45.31 


s 


OJ 


m  ■r  rrj  tv  — • 

•  •  •  •  •  Q 

*NI  ‘30ddans  WOdJ  lN3W3Da3dSia 


u 

Ui 

(ft 

\ 


o 

o 

Ui 

> 

x 

o 


-101- 


v„*.  .-■ 


Fig.  34.  Boundary  Layer  Velocity  Profiles,  Conf.#l  at  50  %  Chord 
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Fig.  40.  Boundary  Layer  Velocity  Profiles,  Conf.#2  at  4.68  %  Chord 
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